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ABSTRACT 

We present high resolution echelle spectra of 7 proximate damped Lyman alpha (PDLA) 
systems whose relative velocity separation from the background quasar is AV < 3000 km s _1 . 
Combining our sample with a further 9 PDLAs from the literature we compare the chemical 
properties of the proximate systems with a control sample of intervening DLAs. The PDLAs 
are usually excluded from statistical studies of absorption selected galaxies and this sample 
constitutes the first systematic study of their chemical and ionization properties. Taken at face 
value, the sample of 16 PDLAs exhibits a wide range of metallicities, ranging from Z ~ 1 /3Z 
down to Z ~ 1/IOOOZq, including the DLA with the lowest N(SiII)/N(HI) yet reported in the 
literature. However, some of these abundances may require ionization corrections. We find 
several pieces of evidence that indicate enhanced ionization and the presence of a hard ion- 
izing spectrum in PDLAs which lead to properties that contrast with the intervening DLAs, 
particularly when the 7V(H i) is low. The abundances of Zn, Si and S in PDLAs with log N(H i) 
> 21, where ionization corrections are minimized, are systematically higher than the interven- 
ing population by a factor of around 3. We also find possible evidence for a higher fraction of 
NV absorbers amongst the PDLAs, although the statistics are still modest. 6/7 of our echelle 
sample show high ionization species (SilV, CIV, OVI or NV) offset by >100 km s -1 from the 
main low ion absorption. We analyse fine-structure transitions of CII* and Sill* to constrain 
the PDLA distance from the QSO. Lower limits range from tens of kpc up to >160 kpc for 
the most stringent limit. We conclude that (at least some) PDLAs do exhibit different char- 
acteristics relative to the intervening population out to 3000 km s _1 (and possibly beyond). 
Nonetheless, the PDLAs appear distinct from lower column density associated systems and 
the inferred QSO-absorber separations mean they are unlikely to be associated with the QSO 
host. No trends with AV are found, although this requires a larger sample with better emis- 
sion redshifts to confirm. We speculate that the PDLAs preferentially sample more massive 
galaxies in more highly clustered regions of the high redshift universe. 
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1 INTRODUCTION 

Absorption systems with 13 < log iV(Hi) < 17 and within ~ 5000 
km s" 1 of a background QSO are often considered to be 'associ- 
ated' with the quasar host galaxy, its immediate environment or its 
ejecta. Evidence for this conclusion comes from time variability 



* Based on observations made with ESO Telescopes at the Paranal Obser- 
vatories under program 080A-0014(A). 



of the absorption (Hamann et al. 1997a; Wise et al. 2004; Lund- 
gren et al. 2007), partial coverage of the emission source (Hamann 
et al. 1997b,c; Fox et al. 2008), elevated incidence of highly ion- 
ized species (e.g. Ganguly et al. 2001; Richards 2001; Fechner & 
Richter 2009), possible dependence of absorber incidence on QSO 
properties such as radio loudness or orientation (M0ller & Jakob- 
sen 1987; Richards et al. 2001 ; Baker et al. 2002; Vestergaard 2003; 
Wild et al. 2008) and solar or super-solar metallicities (e.g. Petit- 
jean, Rauch & Carswell 1994; Hamann et al. 1997a; Srianand & 
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Petitjean 2000; D'Odorico et al. 2004; Fechner & Richter 2009). 
The associated population of absorbers appears to extend out to at 
least 10,000 km s~' (Petitjean et al. 1994; Richards et al. 1999; 
Wild et al. 2008; Tytler et al. 2009). 

Although the z abs ~ z cm Lya forest, CIV, Mgll and Lyman 
limit systems have been well-researched in the last 20 years, the 
same is not true of the proximate damped Lyman alpha systems 
(DLAs). This is partly a historical bias. DLAs are excellent tools 
for characterising the high redshift galaxy population. Knowing 
that associated systems may be connected with the QSO or its 
host galaxy, statistical samples of DLAs have therefore excluded 
these potentially special systems (e.g. Wolfe et al. 1995; Ellison et 
al. 2001; Prochaska, Herbert-Fort & Wolfe 2005; Jorgenson et al. 
2006; Noterdaeme et al. 2009). Furthermore, the chemical abun- 
dances of intervening DLAs can be determined fairly easily be- 
cause their gas is predominantly neutral. Proximity to a QSO could 
invalidate that assumption and lead to erroneous abundance deter- 
minations, unless photoionization modelling is used. However, it is 
for exactly these reasons that the study of proximate DLAs (PD- 
LAs) is interesting. PDLAs provide a potential probe of galaxies 
clustered around (or associated with) QSOs, and even provide in- 
sight into the interplay between ionizing radiation and the ISM. 

One of the difficulties in studying the proximate DLAs is the 
small redshift path for each line of sight, which means that PDLAs 
are very rare. For example, assuming the same number density as 
intervening systems at z = 3, ~ 100 QSOs would need to be sur- 
veyed to find just one PDLA. A second challenge is constructing a 
sample of PDLAs selected with a homogeneous QSO-absorber ve- 
locity offset (AV) criterion. Since the derived emission redshift of 
a QSO depends critically on the emission lines used (e.g. Gaskell 
1982; Tytler & Fan 1992), the value can be uncertain by over 1000 
km s" 1 (~ 3 Mpc proper at z = 3). The first systematic study of PD- 
LAs was conducted with the Complete Optical and Radio Absorp- 
tion Line System (CORALS) Survey data by Ellison et al. (2002) 
who found a factor of 4 excess in the number density of PDLAs 
in their radio-selected sample. Russell, Ellison & Benn (2006) ex- 
tended this work to identify 33 PDLAs in the Sloan Digital Sky Sur- 
vey (SDSS) Data Release 3 (DR3) and confirmed the PDLA excess, 
independent of radio loudness, at 3.5cr significance. In the largest 
PDLA survey conducted so far, Prochaska, Hennawi & Herbert- 
Fort (2008) use the SDSS DR5 to identify 108 PDLAs with AV < 
3000 km s . In addition to its larger size, this work re-calculated 
the QSO redshifts in order to obtain improved estimates of AV. 
Prochaska et al. (2008b) found that PDLAs outnumber the inter- 
vening DLAs at z ~ 3 by a factor of 2, but found no statistically 
significant excess at z < 2.5 or z > 3.5. 

In addition to the PDLA number density, Prochaska et al. 
(2008b) investigated whether their Hi column density distribution 
function differed from the intervening DLAs. They found marginal 
evidence for an excess of high A'(H i) PDLAs. Previously, Tytler 
(1982) had found no excess for the slightly lower N(Hi) Lyman 
limit systems but with a much larger sample Prochaska, O'Meara 
& Worseck (2010) find a deficit. Apart from their incidence and 
A'(H i) distribution, few studies of PDLAs exist to assess whether 
or not they are consistent with the intervening population. In partic- 
ular, PDLAs are rarely targetted for echelle spectroscopy, so little 
is known about their chemical enrichment. Nonetheless, a small 
number of PDLAs have appeared in chemical abundance studies 
(see Table|7]for those that meet our criteria) and they do not appear 
to have strikingly different properties to the intervening DLAs (Lu 
et al. 1996; Akerman et al. 2005; Rix et al. 2007). The range of 
metallicities, a element enhancements and dust depletions are all 



consistent with the range seen amongst intervening DLAs. Rix et 
al. (2007) also found negligible ionization corrections were neces- 
sary for the PDLA in their study, as is usually the case for other 
DLAs. It has therefore usually been concluded that the PDLAs do 
not constitute a distinct population from the intervening systems 
(e.g. M0ller et al. 1998). 

In this work, we present echelle spectra for seven PDLAs 
(Section^, determine improved emission redshifts (Section[3} and 
measure column densities for a range of metal line species (Sec- 
tion |4](. We supplement our sample with other PDLAs critically 
selected from the literature to form a sample of 16 proximate ab- 
sorbers (Section[5j. Ionization properties and metal abundances of 
PDLAs are compared to the intervening systems in Sections [6] and 
[7] Finally, we use measurements of CII* and upper limits to Sill* 
to constrain the distance between the PDLA and the background 
QSO (Section©. 

Unless otherwise stated, we assume Q. M = 0.3, £2a = 0.7 and 
H = 70 km s-' Mpc-'. 



2 SAMPLE, OBSERVATIONS AND DATA REDUCTION 

Proximate DLAs for this study were mostly selected from the 
SDSS DR5 catalogue of Prochaska, Hennawi & Herbert-Fort 
(2008). These are J014049. 18-083942.5, 1014214.74+002324.3, 
Jl 13130+604420, J124020.9+ 145535.6, J160413.9+395121.9 and 
J2321 15.48+142131.5. For simplicity we hereafter refer to these 
QSOs by their SDSS names in the form Ihhmm+ddmm. One 
additional non-SDSS QSO is also in the sample: Q0151+048 
(UM144). SDSS magnitudes are given in Table[T]for all QSOs ex- 
cept Q015 1+048 where we have taken the R-band magnitude from 
the APvQ plate measurements. Our target selection was based on 
relatively bright QSOs with absorbers whose velocity offset was 
found to be AV < 3000 km s" 1 by Prochaska et al. (2008b), and 
with positions accessible during our scheduled observing runs. The 
AV cut-off is somewhat arbitrary and it should be stressed that it is 
unlikely that these velocities correspond to distances inferred from 
the Hubble expansion (d=v/H(z)) if the proximate absorbers are 
linked either to the host galaxy or associated with its gravitational 
potential. In the former case, the velocities may reflect outflows or 
internal motions. In the latter scenario, the AV distribution would 
be indicative of peculiar motions. We provide empirical evidence 
that the value of AV does not correspond to a Hubble distance by 
showing that a number of the PDLAs have negative velocities rela- 
tive to the QSO (Section© and that elemental trends do not corre- 
late with AV (Section©. 

Observations were obtained with either the UV- Visual Echelle 
Spectrograph (UVES) at the Very Large Telescope (VLT) or with 
the High Resolution Echelle Spectrograph (HIRES) on the Keck 
telescope. 

2.1 UVES data acquisition and reduction 

UVES is a dual-arm echelle spectrograph with a grating cross dis- 
perses There are four different cross dispersers available and two 
dichroics (DIC1 and DIC2), which split the incident beam between 
the blue and red arms. The central wavelengths are governed by the 
grating angles and are adjusted to allow the desired spectral cover- 
age to fall on the 3 CCDs (one in the blue arm and two in the red). 

1 www.ast.cam.ac.uk/~apmcat 
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We combined two of UVES' s 'standard' settings (DIC1 346+580 
and DIC2 437+860) in order to obtain spectral coverage over ap- 
proximately 3050 - 9700 A. Integration times are given in Table[T] 
Observations of a ThAr lamp for wavelength calibration followed 
each exposure. The CCDs were binned 2x2 and a 1 arcsecond slit 
was used, resulting in nominal resolutions of R=41,000 and 39,000 
in the blue and red arms, respectively. We observed J0140-0839, 
J0142+0023 and J2321 + 1421 on the nights of 2007 October 2-3 
with UVES in visitor mode. The seeing ranged between 0.6 and 
1 .0 arcsec. 

The spectra were extracted and reduced using the UVES 
pipeline. The extracted spectra from each exposure were converted 
to a common vacuum-heliocentric wavelength scale. Where mul- 
tiple exposures in a given setting were obtained, or where there is 
spectral overlap between settings, the data were combined, weight- 
ing by the inverse of the flux variances. Table [T] lists the typical 
S/N ratios for each target in the final combined spectrum (where 
there is some wavelength overlap between settings). The spectrum 
was normalized by fitting a smooth continuum through unabsorbed 
parts of the data using the Starlink DIPSO packag^. 



2.2 HIRES data acquisition and reduction 

J1240+1455 and J1604+3951 were observed with the Keck/HIRES 
spectrometer (Vogt 1994) using the upgraded CCD mosaic on 2007 
April 27-28. Jl 131+6044 and Q0151+048 were also observed with 
HIRES on 2006 December 25 and 2006 August 17 respectively. 
The first two quasars, which have PDLAs at z abs < 3, were ob- 
served using the blue cross-disperser (HIRESb) with wavelength 
coverage from roughly the atmospheric cutoff to « 6000A. The lat- 
ter two were observed using the red cross-disperser (HIRESr) with 
the kv418 blocking filter giving an effective wavelength coverage 
of 4300-8700A. The echelle and cross disperser angles are given 
in TableQ] along with total exposure times. All of the observations 
were acquired through the C 1 decker which affords a full width at 
half maximum (FWHM) » 6kms~' spectral resolution. The data 
were reduced with the HIRedux pipeline that is available within the 
XIDL software packag^J One significant difference in the analysis 
of the HIRES data, compared to the UVES reductions, is the stage 
at which the the continuum is fitted. Whereas the UVES data are 
normalized as the final reduction step once the echelle orders have 
been combined into one dimension, the XIDL reduction pipeline 
fits a local continuum to each order before the 2-dimensional spec- 
tra are combined into ID. The advantage of this latter approach 
is that the construction of the 1 -dimensional spectrum is simpli- 
fied, e.g. by removing the effects of the blaze. The disadvantage is 
that since damped Lyo- troughs extend over multiple orders, HIRES 
spectra do not lend themselves easily to Lyo- fitting. We discuss this 
further below. 



3 REDSHIFTS 



AV of the QSO. Typical limiting values of AV adopted in the liter- 
ature, when defining samples of associated absorbers, are <3000 to 
<6000 km/s. AV is calculated as: 



AV : 



R 2 - 1 
R 2 + 1 



where c is the speed of light and 



R : 



1 + Zabs 
1 + Zem 



(1) 



(2) 



where z abs and z cm are the absorption and emission redshifts respec- 
tively. It is well known that the determination of QSO emission 
redshift is sensitive to the choice of emission lines (e.g. Gaskell 
1982; Tytler & Fan 1992), where [OIII] XX 4959, 5007 is usually 
considered one of the most reliable indicators of systemic veloc- 
ity. The SDSS composite spectrum of vanden Berk et al. (2001) 
shows that narrow forbidden lines may only be shifted by ± tens 
of km s -1 , but the commonly used hya, CIV and CIII] lines may 
be offset by many hundred km s -1 . The majority of the PDLAs 
used in this study (see Section|2]( are selected from the SDSS DR5 
survey of Prochaska et al. (2008b), and we adopt their re-derived 
emission redshifts. For the QSOs not in the sample of Prochaska et 
al. (2008b) we derive emission redshifts in an identical way. The 
technique first calculates line centres of all significantly detected 
emission lines as described by Hennawi et al. (2006) and applies 
the known systematic shifts of each species as tabulated by Shen et 
al. (2007) and Richards et al. (2002). The high ionization lines such 
as CIV and SilV give the poorest estimate of z em , with a redshift er- 
ror of ~ 700 km s , From the SDSS spectra, our best estimates 
of systemic redshift are derived from the Mgll emission line which 
yields z em accurate to ~ 270 km s , based on the distribution of 
MgII-[OIII] offsets given in Richards et al. (2002). Emission red- 
shift errors in Table [2] are based on the distributions of line shifts 
given in Table 7 of Shen et al. (2007). For future studies, it would 
be desirable to obtain infra-red (IR) spectra of these targets to ob- 
tain redshifts directly from the Balmer and forbidden [OIII] lines 
(e.g. Rix et al. 2007), the latter of which yields accuracies of less 
than 50 km s . Recently Hennawi et al. (2009) have obtained an 
IR spectrum for one of our sightlines: 11240+1455 and for this ob- 
ject we adopt the emission redshift determined from the [OIII] line 
(Zcm = 3.1092). This value is larger by 160 km s _I than the value 
derived from the SDSS spectrum (z em = 3.1070) and is now larger 
than the absorption redshift. 

Having made our best estimate of the emission redshifts of the 
QSOs in our sample, a value of the absorption redshift must also be 
obtained. The velocity spread of metal species in DLAs can extend 
over several hundred km s _I (e.g. Prochaska et al. 2008a) and the 
redshift of the best fitting Lyo- line can be offset from the metals. In 
Table|2]we give the absorption redshifts of our sample determined 
from both the Lyo- line and the strongest metal component (z a b s ra 
and z absZ respectively). It can be seen that the resulting difference 
in AV is much smaller than the typical error associated with the 
emission redshift determination. Since the Lyo- redshift is the value 
most easily recovered from the literature, we adopt that as our fidu- 
cial measure of z abs . 



Accurate emission and absorption redshifts are required in or- 
der to define a sample of PDLAs lying within a given velocity range 

2 http://star-www.rl.ac.uk/staiydocs/sun50.htx/sun50.html 

3 http://www. ucolick.org/~xavier/IDL 



4 COLUMN DENSITY MEASUREMENTS 

For the 3 targets with UVES data, Hi column densities were de- 
rived by fitting damped profiles to the Lyo- absorption in the nor- 
malized spectra using the DIPSO software package. Continuum 
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QSO rmag z cm Zabs log /V(H i) (cm 2 ) Instrument Instrument Exposure S/N 

(SDSS) (UVES/SDSS) set-up time (s) pix" 1 



J0140-0839 


17.7 


3.7156 


3.6960 


20.75±0.15 


UVES 


346+580 


6000 


30- 


50 














437+860 


6000 






J0142+0023 


18.3 


3.3734 


3.34765 


20.38±0.05 


UVES 


346+580 


6000 


30- 


50 














437+860 


10200 






Q0151+048 


17.5 


1.9225 


1.9342 


20.34±0.02 


HIRES 


0/1.250 


9000 


15- 


20 


J1131+6044 


17.7 


2.9069 


2.8754 


20.50±0.15 


HIRES 


0/1.680 


7200 


15- 


20 


J1240+1455 


18.9 


3.1092 


3.1078 


21.3±0.2 


HIRES 


-0.110/0.557 


5400 


5- 


8 


J 1604+3951 


18.1 


3.1542 


3.1633 


21.75±0.2 


HIRES 


-0.110/0.557 


10300 


10- 


15 


J2321+1421 


18.3 


2.5539 


2.5731 


20.70±0.05 


UVES 


346+580 


3000 


10- 


30 














437+860 


11700 







Table 1. Target list and observing journal. The instrument set-up for HIRES refers to the echelle/cross-disperser angles. The S/N is given as a representative 
value of the final combined spectra of both settings. 



QSO 


log N(HI) 


-em 


ZabsHI 


ZabsZ 


AV H i 


AV Z 




(cm" 2 ) 








(km s _1 ) 


(kms^ 1 ) 


J0140-0839 


20.75±0.15 


3.7156±0.012 


3.6960 


3.69660 


1250 


1211 


J0142+0023 


20.38±0.05 


3.3734±0.008 


3.34765 


3.34768 


1772 


1769 


Q0151+048 


20.34±0.02 


1.9225±0.003 


1.9342 


1.93429 


-1199 


-1208 


Jl 131+6044 


20.50±0.15 


2.9069±0.009 


2.8754 


2.87562 


2424 


2412 


J1240+1455 


21.3±0.2 


3.1092±0.005 


3.1078 


3.10803 


102 


85 


J 1604+3951 


21.75±0.2 


3.1542±0.007 


3.1633 


3.16711 


-656 


-930 


J2321+1421 


20.70±0.05 


2.5539±0.006 


2.5731 


2.57312 


-1616 


-1618 



Table 2. PDLA redshifts and relative velocities 



J0140-0839 z ab .=3.6760, ".og K(HI) - 20.75 
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Figure 1. Fits to the Lyo- profiles for the three PDLAs in our sample ob- 
served with UVES. Column densities and uncertainties are given in Table 

m 



adjustment and H i fitting were carried out iteratively, due to prox- 
imity of the absorption line to the strong Lycr emission of the QSO 
itself. The Hi column densities derived from the UVES data are 
all in good agreement with those derived from the SDSS spectra 
published in Prochaska et al. (2008b). Figure [T] shows the Lyo- fits 
to the PDLAs. For 2/3 of the UVES targets we were able to check 
our Lyo- fits against an independent fit of Ly/? (for JO 142+0023 Ly/J 
is blended with strong Lyo- at a lower redshift). In both cases the 
Lya and LyyS column densities agree to within the errors quoted in 
Tables[j]and[2] 

For the 4 HIRES targets where DLA fitting is hampered by 
the order-by-order normalization in XIDL, we adopt the HI column 
density derived from the SDSS spectrum (Prochaska et al. 2008b). 
As described above, the fits to the UVES spectra demonstrate that 
the SDSS values agree with the higher resolution spectra fits to 
within the errors. For Q015 1+048, the H i column density was mea- 
sured from an X-shooter spectrum by Zafar et al. (in preparation), 
which agrees with the value of M0ller et al. (1998) to within 0.02 
dex. 

Metal column densities were mostly derived by fitting Voigt 
profiles with the VPFIT packagfl The fitting of Voigt profiles is 
particularly useful in distinguishing contamination from blended 
absorption features and simultaneously fitting multiple transitions 
of a given species. In a few cases (noted below), usually when 
the line was too weak, or the S/N too low for a convincing Voigt 
profile decomposition, the apparent optical depth method (AODM, 
e.g. Savage & Sembach 1991) was applied. In general, the kine- 
matic absorption model parameterized by redshifts and fc-values 
(Doppler widths) for multiple components was derived indepen- 
dently for each species. Again, exceptions to this are noted below, 



4 http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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usually in cases where the lines are weak or blending is suspected. 
Tables [3] and [4] list the adopted column density for low and high 
ionization species respectively. 

The calculation of upper limits depends on an assumption of 
the number of pixels over which a line should be detected. For a 
line of given FWHM in A, the observed frame n sigma equivalent 
width (EW) detection limit is given by 

n x FWHM 

EW = (3) 

S/N V ' 

where S/N is the signal-to-noise per pixel. In order to determine 
the appropriate value of the FWHM, we take the fo-value of the 
strongest component of a detected species for a given line of sight 
(this is usually a transition of Sill or Fell) and apply the correction 
FWHM=£> x 2 V'i2, before converting from km s~ 1 into A. A single 
value of the FWHM is used for all low ions, CII* and A1III limits in 
a given sightline. Typical i-values are 4-8 km s -1 . If the strongest 
component is unresolved, a FWHM corresponding to the instru- 
mental resolution is adopted, i.e. the FWHM that appears in Eqn 
[3] is the maximum of either the strongest component or the instru- 
mental resolution. Had we assumed that the lines were unresolved 
rather than adopting true line widths, the upper limits would have 
been under-estimated by 0.2-0.3 dex. For CIV, SilV and NV we 
assume a broader profile and adopt a single value of b= 10 km s . 
Upper limits in this paper are quoted at 3<x significance. 

For species where only saturated lines are detected, lower limit 
column densities have been derived by applying the AODM to the 
observed line with the lowest oscillator strength (/ value). Column 
densities are converted to abundances by adopting the solar scale 
of Asplund et al. (2005), with the exception of argon for which 
we adopt the improved value of Asplund et al. (2009). Meteoritic 
values are adopted with the exception of C, N and O where we 
take the photospheric values (and Ar where several indirect meth- 
ods are used, see Asplund et al 2009 for details). We note that there 
are significant differences in the solar abundances of some of these 
volatile elements relative to other commonly used reference scales, 
e.g. Holweger (2001). All literature abundances used for compar- 
ison in the present work have therefore been re-calculated on the 
solar scale of Asplund et al. (2005). Following the standard pro- 
cedure in DLA abundances, we assume that the singly ionized 
species dominate the total column densities (i.e. for element X, 
N(XII)=N(X)) and we do not apply ionization corrections. Even 
though this assumption may not hold in some of the PDLAs in 
our sample (we discuss ionization in Section |6j, it is nonetheless 
of interest to apply this common practice to the PDLAs to assess 
its impact. Exceptions to this approximation are three-fold for the 
elements reported from our data. Charge exchange reactions (and 
ionization potentials < 13.6 eV) keep the majority of nitrogen and 
oxygen in the neutral state and for aluminium both A1II and A1III 
can contribute significantly. Since we do not calculate ionization 
corrections on a system-by-system basis, aluminium is omitted in 
our final abundance listing, see Table[6] 

Abundances are quoted relative to solar using the usual nota- 
tion 



[XI H] = log 



N(X) 



N(H) 



log 



MX) 



N(H) 



(4) 



4.1 J0140-0839, log N(H i) = 20.75 

The PDLA towards J0140-0839 has a simple velocity structure 
with two main components in the low ionization species at v =0 
and +50 km s (see Figure [2). Three unsaturated transitions of 



Sill are detected (A =1260, 1304, 1526 A), although the latter of 
these exhibits some blending. The column density is derived by 
fitting the two bluer transitions simultaneously in VPFIT. Arl A 
1048 is blended but the redder transition (Arl A 1066) of the dou- 
blet coincides with a clean region of the spectrum and is used to 
determine an upper limit to the argon column density. Despite the 
fairly high S/N (30 - 50), A1II is barely detected and its column 
density is derived using both VPFIT and the AODM. In VPFIT we 
adopt the fo-values and redshifts of the two strongest components 
of the Sill fit to derive log N(A1II)= 1 1 .78±0.04. The AODM gives 
log N(A1II)=11.85±0.14, in agreement with the fit, within the er- 
rors. We adopt the average value of these two methods. No high 
ionization species were detected for this PDLA. CIV is usually de- 
tected even in the low metallicity Lyce forest (e.g. Cowie et al. 1995; 
Songaila 1997; Ellison et al. 1999, 2000), and at column densities 
1-3 orders of magnitude above our detection limit of N(CIV) < 
12.18. However, this PDLA is extremely metal-poor; both its Si 
and Fe abundances are the lowest in the current published literature 
(of which we are aware). Fox et al. (2007) have shown that there is a 
broad correlation between N(CIV) and metallicity, consistent with 
the non-detection of CIV in this PDLA. As we discuss in Section 
[6]there may be corrections due to ionization effects, (and addition- 
ally some dust depletion, at least in the case of Fe). However, the 
OI is little affected by dust and ionization and the OI A 1302 line is 
unsaturated in this PDLA, from which we determine [0/H]=-2.72. 
This supports an intrinsically low metallicity for this system. In- 
deed, the [O/H] of this PDLA is similar to the oxygen abundance 
of the IGM (Simcoe et al. 2004; Aguirre et al. 2008) who find [O/H] 
2.5 to -3 depending on the local over-density, although this de- 
pends on the model for ionization corrections and shape of the UV 
background. Previously, a metallicity 'floor' has been suggested for 
DLAs (e.g. Prochaska et al. 2003a) and this system offers a rare 
chance to measure abundances at extremely low levels of enrich- 
ment. For example, Akerman et al. (2005) and Pettini et al. (2008) 
discuss how carbon abundances at low metallicity might provide 
clues into the yields of early stellar populations. The unsaturated 
CII A 1334 yields [C/O] =-0.29 supporting the idea of an upturn in 
C/O ratios at low (O/H) (Pettini et al. 2008). The non-detection of 
N also sets a very low upper limit to N/O: [N/O] < -1.43. 



4.2 J0142+0023, log N(R i) = 20.38 

Selected metal lines are shown in Figure [3] Two lines of Sill (A 
1304 and A 1526) are detected, although they are approaching satu- 
ration. N(SiII) was obtained separately for the two transitions with 
VPFIT and found to be in excellent agreement. A1II A 1670 fol- 
lows the velocity structure of the other low ions. There is weak 
absorption at A1III A 1 862 at a similar velocity, although gas in a 
different ionization state need not necessarily trace the neutral gas. 
A comparison with A1III A 1854 could confirm the detection, but 
unfortunately this line is blended, so the detection of the A 1862 
line remains uncertain. We therefore calculate N(AIIII) from A1III 
/U862 and take this as an upper limit. SII A 1259, although weak, is 
statistically a 5cr detection. Absorption is detected from both SilV 
and CIV. Absorption is present in a similar velocity range as the low 
ions, but with additional components present at negative velocites 
from -150 to km s" 1 , There is a significantly detected feature at 
the expected velocity of Arl A 1048, but additional absorption to the 
red alerts us to the possibility of contamination. The identification 
of this line as Arl A 1048 could be confirmed by the simultaneous 
detection of Arl A 1066, but this line is blended with Lyo-. Assum- 
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Figure 2. Selected metal line transitions in the PDLA towards J0140-0839. The lower solid (red) line shows the error array. When fits have been derived using 
VPFIT, those fits are overlaid in green, b- values for this absorber range from 2-8 km s . Velocities are plotted relative to z a b s =3.6966. 



ing the v = feature is Arl A 1048, we derive a column density 
from the AODM and quote this as a conservative upper limit. 



4.3 Q0151+048, log N(B i) = 20.34 

Selected metal lines for the PDLA towards Q015 1+048 are shown 
in Figures[4]and[5] As described by M0ller et al (1998) and Fynbo 
et al. (1999, 2000), Q0151+048 has a close, but fainter, compan- 
ion, Q0151+048B, which is separated by 3.27 arcsec (~ 27.5 
kpc). M0ller et al. (1998) give the redshift of Q0151+048 to be 
Zem = 1-922 ± 0.003 and a slightly higher value for Q0151+048B 
of Zem = 1-937 ± 0.005. As shown by Zafar et al. (in prepara- 
tion) Q0151+048B shows no DLA absorption, despite its rela- 
tively small transverse separation. Using the same technique de- 
scribed in Section [3] we re-evaluate the redshift of Q0151+048 
and Q0151+048B from X-shooter spectra kindly provided by Jo- 
han Fynbo and Tayyaba Zafar. We determine Zcm = 1.9225±0.003 
for Q0151+048 and z era = 1.9237±0.003 for Q0151+048B, where 
both redshifts are based on the Mg II line. However, we caution that 
the spectrum of Q0151+048B is affected by a lot of structure in the 
continuum shape which may compromise the accuracy of the red- 
shift (the quoted error only accounts for known systematic offsets 
of emission lines from the systemic value). Interestingly, the PDLA 
studied by Rix et al. (2007) towards Q2343-415 also has a nearby 
QSO companion (Q2343+125) separated by 680 kpc and -100 
km s . 



N(FeII) is derived from a simultaneous VPFIT of Fell A 1608 
and Fell A 1 144. Similarly, N(SiII) is determined by simultaneously 
fitting the lines at A 1304 and 1526 A. The NI column density was 
determined from a simultaneous fit of NI A 1 199 and NI A 1200.2. 
CII* is detected, but its low optical depth precludes an accurate 
Voigt profile fit; its column density is derived from the AODM. 
By combining N(CII*) with A'(Hi) it is possible to determine the 
cooling rate, l c , which Wolfe et al. (2008) have recently shown to 
be bi-modal. We calculate a cooling rate log l c = -26.86 which 
qualifies this PDLA as a 'high cool' (i.e. high cooling rate) sys- 
tem. Such systems typically have higher metallicities, dust-to-gas 
ratios and velocity spreads than the low cooling rate DLAs. How- 
ever, the PDLA towards Q0151+048 has a relatively narrow veloc- 
ity structure and apparently low metallicity ([S/H]< -2.03). Either 
this PDLA is rather unusual amongst the 'high cool' population, 
or the column densities may be affected by ionization. However, it 
should be noted that it lies towards the lower end of the l c distribu- 
tion for high cool systems. The low N(H i) of this system make the 
latter possibility a likely explanation and we present more evidence 
for ionization in low A'(Hi) PDLAs in Section|6] 

Both CIV and SilV are very strong and may be partly satu- 
rated. The velocity structure is simple and coincides with that of 
the low ions, see Figure [4] which is usually not the case for inter- 
vening DLAs. AODM and VPFIT give consistent lower limits for 
CIV. N(SiIV) is derived from the A 1402 line due to mild satura- 
tion in the bluer member of the doublet. In Section [6] we discuss 
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Figure 3. Selected metal line transitions in the PDLA towards JO 142+0023. The lower solid (red) line shows the error array. When fits have been derived using 
VPFIT, those fits are overlaid in green, b- values for this absorber range from 3-9 km s for the low ionization species and 7-25 km s -1 for the higher 
ionization species. Velocities are plotted relative to z a bs= 3.3477. 



the relative contributions of Sill and SilV in this PDLA and the 
implications of coincident velocity structure. 



Although no NV is found at the redshift of the main absorber, 

there is NV offset by 825 km s -1 . This velocity corresponds to a 

redshift that is well offset into the Lyo- wing, see Figure[5] and does 
not appear to be associated with high column density HI. The NV 
is well-fit in VPFIT (Figure [5] so there is no evidence for partial 
coverage in the current data. The NV is accompanied by weak CIV, 
but no SilV (OVI is not covered by the spectra). The column densi- 
ties of this component are not included in the main table (Table [6} 
as this absorption does not appear to be associated with the same 
PDLA system. From VPFIT we determine N(NV)=13.12±0.03 
from the U 1238, 1242 doublet and N(CIV)=12.89±0.02. This 
offset absorption is qualitatively similar to the intervening systems 
studied by Schaye et al. (2007), who describe a population of ab- 
sorbers with significant column densities of CIV and NV, but with 
little HI. Weidinger et al. (2005) also find orphaned NV near a prox- 
imate Lyman limits system with even larger velocity offset. Fechner 
& Richter (2009) give the line of sight number density of interven- 
ing NV absorbers to be dN/dz=1.16. For a Az = ±0.01 = 0.02 
(which corresponds to ± ~ 1000 km s _I at z = 2) there is a 0.2% 
random chance that an intervening NV system lies this close to a 
PDLA. We discuss the offset NV further in Section|9] 



4.4 J1131+6044, log N(H i) = 20.50 

Selected metal line transitions for the PDLA towards Jl 131+6044 
are shown in Figure(6] The Sill column density is determined solely 
from the A 1526 line, due to saturation of Sill A 1260, blending of 
the Sill A 1304 line and no coverage of Sill A 1808. Similarly, the 
only detected Fell line is at /lo=1608 A. Fixing the velocity struc- 
ture to match that of the Sill fit produces an excellent fit. Of the 
NI triplet, only the bluest transition at /l =1199 A is unblended. 
We determine N(NI) using both the fixed component model de- 
rived from Sill and using a free fit. The former produces a rel- 
atively poor fit (but does not indicate Lyo- blending) to the data 
and yields N(NI)=13.87±0.09. The latter yields a slightly lower 
N(NI)=13.76±0.08. We adopt an intermediate value of 13.8+0.15. 
There is a weak (5 a) feature at the expected position of Arl 1048 
whose column density can be determined from the AODM. Ideally, 
detection of Arl A is required to confirm that this is not contami- 
nation by weak Lya. However, since Arl A 1066 is not detected to 
confirm the identification, we take the conservative approach and 
use the AODM-derived column density as a limit. 

Although this PDLA has a relatively high AV (~ 2400 km s~') 
it still demonstrates some interesting properties, most notably in 
its high ionization lines. At v ~ the CIV and SilV absorption is 
relatively weak. The SilV AA 1393, 1402 doublet is fitted simult- 
neously. The CIV /U548 line is blended with extended negative ve- 
locity gas (see below and Figure[7} so N(CIV) is determined solely 
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Figure 4. Selected metal line transitions in the PDLA towards Q015 1+048. The lower solid (red) line shows the error array. When fits have been derived using 
VPFIT, those fits are overlaid in green, b- values for this absorber range from 2-8 km s for the low ionization species and 5-15 km s -1 for the higher 
ionization species. Velocities are plotted relative to Zabs= 1.9342. 



from CIV A 1550. There is no NV at v = but we are not able 
to determine whether or not OVI is present due to the blending 
in at least one of the doublet components. However, the kinematic 
similarity between absorption at the expected wavelength of OVI/1 
1031 and CIV A 1550 is suggestive (Figure[7J. We do find highly 
offset absorption of strong CIV and OVI (but no NV) at v - — 750 
km s -1 . Since we can not estimate the total OVI column density due 
to blending, no value is given in Table H The CIV extend^ from 
approximately -500 to - 1 150 km s . We can not confirm the full 
extent of the OVI in velocity space due to various hya blends, see 
Figure [7] The metal column densities in this extended component 
are not included in the values in Table [4] 



4.5 J1240+1455, log N(H i) = 21.3 

Selected metal lines are shown for the PDLA towards J1240+1455 
in Figured] This PDLA was studied by Hennawi et al. (2009) using 
the original SDSS spectrum. Like several other PDLAs (M0ller & 
Warren 1993; M0ller et al. 1998; Leibundgut & Robertson 1999; 
Ellison et al. 2002) this absorber has Lya emission superimposed 



5 There are two weakly detected features in this approximate velocity range 
that could be NV /11238. However, NV A 1242 is not detected at these 



extended velocites. 
present data. 



so these identifications can not be confirmed with the 



in the DLA trough. After extensive modelling, Hennawi et. (2009) 
concluded that the Lya emission, is likely to be associated with the 
QSO, and not the PDLA itself. 



This is a relatively faint QSO and the HIRES spectrum ex- 
hibits a low S/N. Nonetheless, we are able to detect a number of 
metal lines due to the high V(Hi) (see Figure [8}. For example, 
all three lines in the SII AAA 1250, 1253, 1259 triplet are detected 
and are fitted simultaneously. We check for saturation by addition- 
ally calculating the AODM column densities, which yield the same 
column density in each triplet. Due to the low S/N, we determine 
N(FeII, Sill) in three different ways: using a fixed component VP- 
FIT model based on the component structure of SII, a model where 
the component structure can vary and finally the AODM. All give 
very consistent answers, so we adopt that of the fixed component 
VPFIT model. Since the optical depths in the line centres of the 
Sill and Fell lines are less than in the demonstrably unsaturated SII 
lines, saturation does not seem to be an issue. The Fe and Cr abun- 
dances are both low in this system (e.g. relative to Zn), indicating 
that dust depletion may be significant. Strong, broad NV absorp- 
tion is detected, offset by ~ -120 km s from the strongest low 
ion component. Only the weaker NV A 1242 line is fitted, although 
even this may be partly affected by saturation. We therefore conser- 
vatively quote the fit as a lower limit in Table [4] CIV and SilV are 
both heavily saturated and CII* is blended. 
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Figure 8. Selected metal line transitions in the PDLA towards J1240+1455. The lower solid (red) line shows the error array. When fits have been derived using 
VPFIT, those fits are overlaid in green, fe-values for this absorber range from 7-18 km s -1 for the low ionization species and 6-30 km s for the higher 
ionization species. Velocities are plotted relative to z a bs= 3.1078. 



4.6 J1604+3951, log N(H i) = 21.75 

The PDLA towards J 1604+3951 exhibits two main components in 
the low ionization gas separated by about 250 km s , see Figure[9] 
However, the strongest transitions, such as CII A 1334 and A1II A 
1670 show almost continuous absorption over 600 km s _1 . Singly 
and multiply ionized species exhibit quite similar velocity struc- 
ture, e.g. Sill and SilV (Figure [TOt. In the alpha elements (Si and 
S) the redder component is stronger, but in Fe-peak (Fe and Ni) 
the blue component is stronger. The absorption structure of Znll 
is more akin to Si and S than Fe and Ni, even though it is often 
considered to track the Fe-peak. We return to this point below. 

SII A 1250 and A 1253 are fitted simultaneously and we apply 
the same model for Sill A 1808. Fell /11608 is saturated in several 
components and yields alower limit of N(FeII)>15.15. Fortunately, 
Fell /1161 1 is detected, albeit with low optical depth. We try both a 
fixed (tied to the SII structure) and free velocity model; both give 
a consistent answer of N(FeII)= 15.35+0.03. The AODM gives a 
slightly higher value of N(FeII)= 15.45+0.2. We adopt an interme- 
diate value of N(FeII)= 15.4+0. 15. CIFM1335 is saturated at v ~ 
km s -1 and the negative velocity component is blended with CII 
A 1334. We therefore quote a lower limit that is very conservative 
since it does not (cannot) account for the negative velocity gas. Sat- 
uration and blending mean that the lower limit for CII A 1334 is 
similarly conservative. As for Fell, we attempt three different ap- 
proaches for Zn, all give consistent values, so an average is taken. 



Arl A 1048 is blended, as are the weaker components of Arl A 1066. 
However, the strongest velocity component of Arl A 1066 appears 
clean. Whilst it is therefore not possible to report a total N(Arl), 
we fit the Arl A 1066 line with the same velocity structure as Sill 
to derive N(ArI)= 14.45 for the v = complex. The Sill column 
density for these components is 15.96, yielding [Ar/Si]=-0.40. 

Figure \10\ shows absorption from the more highly ionized 
species. CIV and SilV are particularly strong. Due to the wide 
velocity structure, CIV is self-blended. The limit for N(CIV) is 
therefore derived from the blue (negative velocity) component of 
CIV /11548. NV is detected, but the poor continuum determina- 
tion makes fitting Voigt profiles challenging. A simultaneous fit 
to NV ,11238 plus NV A 1242 from -100 to 250 km s" 1 yields 
N(NV)= 14. 14+0.02 (the error does not include errors in the con- 
tinuum). There appears to be further NV at -400 km s" 1 , where 
CIV and SilV are also seen, but not the low ionization species. Fits 
to this bluer component are hampered by the poor continuum place- 
ment and possible blending from other lines and not accounted for 
in the column density quoted in Table|4] 

Considering the relative abundances produces some clues to 
the difference in the velocity structure of the low ions. The col- 
umn densities for the two main velocity components at -240 and 
km s -1 are given in Table[5] Also tabulated are the relative abun- 
dances of some key elemental ratios. Both components have [S/Si] 
ratios consistent with the solar value, as expected from studies 
of Galactic stars (Chen et al. 2002). Similarly, the [S/Zn] ratio is 



10 Ellison et al. 



1.5 
1 

0.5 



0.5 



Znll 2026, 




1 

x 

E 0.5 

^ 

-400-200 -400-200 -400-200 -B 006001 04)2000 200 
Velocity (krn/s) Velocity (km/s) Velocity (km/s) Velocity (km/s) 

Figure 9. Low ionization metal line transitions in the PDLA towards J 1 604+395 1 . The lower solid (red) line shows the error array. When fits have been derived 
using VPFIT, those fits are overlaid in green, fo-values for this absorber range from 3-12 km s" 1 for the low ionization species. Velocities are plotted relative 
to z a bs= 3.1670. Note the different velocity scale of the last column. The CII* line is blended to the blue with CII. 



consistent with mild a element enhancement in both components 
(where S and Zn are both undepleted and therefore a useful com- 
bination for this assessment, e.g. Nissen et al. 2004). However, 
whereas the blue component has a solar [Zn/Fe] ratio, the gas at 
v ~ km s -1 has a very high value: [Zn/Fe]=0.75. These results are 
strongly indicative of depletion patterns varying within the galaxy. 

There is a possible detection of SiII*/i 1264. Although ex- 
cited transitions of Si and Fe have been previously reported in 
DLAs associated with GRBs (e.g. Prochaska, Chen & Bloom 2006; 
Vreeswijk et al. 2007), this is the first possible detection in a QSO 
(P)DLA. The fit gives N(Sin*)=12.6±0.1. Although this column 
density is 16 times larger than the 3<x upper limit derived for the 
absorber PKS 1443+27 (Howk et al. 2005), the iV(Hi) of the PDLA 
is also almost 10 times higher. Unfortunately, the saturation of the 
CII* transition precludes an estimate of the warm neutral medium 
fraction in this absorber. 



4.7 J2321+1421, log N(H i) = 20.70 

The PDLA towards J2321 + 1421 has a relatively simple velocity 
structure with most of the absorption centred within ±40 km s _1 of 
z a bs = 2.5731, see FigureQT] N(FeII) is determined from a simulta- 
neous fit to Fell A 1608 and Fell A 2374, both of which are unsatu- 
rated. Sill A 1808 is a maginal detection, whereas Sill A 1304 and 
Sill A 1526 are mildly saturated. Nonetheless, combining the 3 tran- 
sitions we are able to constrain the total Sill column density. The 



region of spectrum around the strongest CiTI line at /lo=2056 A is 
relatively noisy, so we derive an upper limit from the non-detection 
of CrII A 2062. There is minor blending of NI with Lyo- forest lines 
but simultaneous fitting of the triplet gives a consistent fit. There 
is weak (6 a) absorption at approximately the expected position of 
A1III A 1854; N(A1III)=12.2 is determined from the AODM. How- 
ever, this is inconsistent with the non-detection of A1III 1862 for 
which we derive a 3 it limit of N(A1III)<11.9. Given the weakness 
of the 1854 A line and the possibility of contamination we give 
12.2 as a conservative upper limit. CIV and SilV are detected, but 
offset by -100 km s -1 from the strongest low ion absorption, but 
coincident with much weaker column density components e.g. in 
A1II and Fell (see Figure ITTb . The spectral region around the Arl 
doublet is fairly noisy. Arl A 1066 is not detected, but there is a 5 
cr feature at the expected velocity of Arl A 1048. However, the S/N 
at this wavelength is only around 2 per pixel, so we consider this 
unreliable and thterefore conservatively adopt the upper limit from 
Arl A 1066. 
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Figure 10. High ionization metal line transitions in the PDLA towards J 1604+3951. The lower solid (red) line shows the error array. When fits have been 
derived using VPFIT, those fits are overlaid in green, b- values for this absorber range from 3-25 km s for the high ionization species. Velocities are plotted 
relative to z a bs= 3.1670. The red component of CIV 1548 is blended with the blue component of CIV 1550. The upper row shows some of the low ions 
(repeated from Figure[9) for comparison. 
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log N(A1II) 
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J0140-0839 


<12.73 


14.13±0.08 


13.51±0.09 


<13.33 


14.69±0.01 


blend 


<12.39 


<12.38 


11.82±0.04 


<12.82 


JO 142+0023 


13.7+0.1 




14.15±0.03 


13.28±0.06 


>15.01 


<11.50 


<12.17 


blend 


12.73±0.01 


<12.57 


Q0151+048 


13.70±0.01 


> 14.43 


14.01±0.05 


< 13.47 


>14.84 


<11.81 


< 12.45 


13.06±0.05 


12.57±0.05 




Jl 131+6044 


13.76±0.03 


>14.55 


14.49±0.13 


< 13.29 


>14.82 






13.8+0.15 




<12.52 


J 1240+ 1455 


14.60±0.03 


sat/blend 


15.93±0.03 


15.56±0.02 


>15.24 


12.90±0.07 


<13.02 


blend 


> 13.56 




J 1604+3951 


15.40±0.15 


>15.28 


16.09±0.02 


15.70±0.02 


blend 


13.0±0.1 






>14.00 


yes/blend 


J232 1+1421 


14.18+0.03 


> 14.68 


14.45±0.04 


<13.60 


>15.10 


<1 1.84 


<12.57 


13.64±0.03 


12.99±0.02 


<13.33 



Table 3. Column densities (in cm 2 ) of low ionization species. 



QSO 



logN(CIV) logN(SilV) logN(NV) log N(OVI) log N(AIIII) N(CII*) 



J0140-0839 


<12.18 


<12.20 


<12.96 


blend 


< 11.52 


< 12.41 


J0142+0023 


14.25±0.01 


13.73±0.01 


<12.29 


blend 


< 12.4 




Q0151+048 


>14.50 


13.75±0.01 


<12.66 




12.3±0.1 


13.0±0.2 


J1131+6044 


13.85±0.05 


13.33±0.02 


<12.68 


yes/blend 




< 12.51 


J1240+1455 


>15.13 


>14.31 


>14.86 




blend 


blend 


J 1604+3951 


>15.05 


>14.74 


14.14±0.02 


blend 


13.44±0.02 


>14.30 


J2321+1421 


13.81±0.05 


13.37±0.01 


<12.62 


blend 


<12.2 


< 12.55 



Table 4. Column densities (in cm ) of high ionization and excited species 
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Figure 11. Selected metal line transitions in the PDLA towards J2321 + 1421. The lower solid (red) line shows the error array. When fits have been derived 
using VPFIT, those fits are overlaid in green, fa-values for this absorber range from 2-12 km s~' for the low ionization species and 5-25 km s -1 for the 
higher ionization species. Velocities are plotted relative to z a bs= 2.5731. The vertical ticks in the lower left panel show the strongest components in each of the 
NI triplet lines. 
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> 


-2.86 








-1.17 




-1.21 


-1.80 




1.36 










J2321+1421 


20.70±0.05 


> 


-2.41 


-2.84 


> 


-2.26 


-1.76 


< 


-2.26 


-1.97 


< 


-1.47 


< 


-1.76 


< 


-1.77 



Table 6. Abundances relative to the solar scale of Asplund et al. (2005), except for argon for which we use the value given in Asplund et al (2009). The 
approximation of N(X)= N(XII) or N(XI) has been assumed. 



5 ADDITIONAL PROXIMATE AND INTERVENING 5.1 PDLAs in the literature 

DLAS FROM THE LITERATURE 



Although the study presented here represents the first systematic 
study of a sample of PDLAs, a small number of proximate systems 
have been included in literature studies of intervening absorbers. In 
this section, we search the literature for PDLAs that meet our AV 
selection criterion and with measured metal column densities that 
can be used to enlarge our sample. We also describe the compila- 
tion of a comparison sample of intervening DLAs with which the 
PDLAs can be compared. 



We searched for additional PDLAs with known abundances by 
trawling the catalogue of Dessauges-Zavadsky et al. (in prep.). 
Emission redshifts were taken from the references given in 
Dessauges-Zavadsky et al. (in prep); if no emission redshift was 
present in the referenced paper, we used SIMBAD. Due to the nu- 
merous different techniques and inherent uncertainties in emission 
redshift determination, a first cut was made for literature DLAs 
with AV < 5000 km s , of which there are 15 in the Dessauges- 
Zavadsky catalogue. We obtained literature spectra for 14/15 of 
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Figure 5. Highly offset NV and CIV in the PDLA towards Q015 1+048. The 
lower solid (red) line shows the error array. When fits have been derived 
using VPFIT, those fits are overlaid in green. Velocities are plotted relative 
toz abs = 1.9342. 





Blue cpt. 

v ~ -240km s _1 


Red cpt. 

v ~ km s 


N(ZnII) 


12.43 


12.88 


N(FeII) 


15.27 


14.97 


N(SiII) 


15.50 


15.96 


N(SII) 


15.16 


15.55 


[Zn/Fe] 


0.00 


+0.75 


[S/Si] 


0.01 


-0.06 


[S/Zn] 


0.18 


0.12 



Table 5. Column densities (cm 2 ) and abundance ratios in two components 
in PDLA J 1604+3951 

these targets that were suitable for a re-assessmen using the 

technique described in Section[3] In this way, all of the PDLA emis- 
sion redshifts and hence the AV values are computed in a consistent 
way. In Table [7] we list the PDLA candidates selected from the lit- 
erature together with updated emission redshifts and re-computed 
AV where available. There is no obvious trend of AV ncw - AV old 



6 We were unable to locate a suitable spectrum for the re-measurement of 
the emission redshift of Q0201+365. 
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Figure 6. Selected metal line transitions in the PDLA towards J 1 1 3 1 +6044. 
The lower solid (red) line shows the error array, fo-values for this absorber 
range from 2-8 km s~' for the low ionization species and 2-10 km s~' 
for the higher ionization species. When fits have been derived using VP- 
FIT, those fits are overlaid in green. Velocities are plotted relative to z a b s = 
2.87562. 



with redshift and positive and negative changes are present in ap- 
proximately equal number. Only one QSO's redshift changes by 
> 1000 km s" 1 , name Q0425-5214 (CTS 436) whose redshift is 
obtained from the Lyo- line (Maza et al. 1995). Not only is this 
line a poor indicator of systemic redshift, but the spectrum is of 
very low disperison (30A resolution) are reported as being accurate 
to only ±0.02. We discard absorbers with AV > 3000 km s -1 . To 
complement the abundances determined for our new data sample 
(Section[4j column densities for literature PDLAs with AV < 3000 
km s -1 are taken from Dessauges-Zavadsky et al. (in prep.) and 
combined with our new echelle sample. The final PDLA sample 
therefore contains 7+9 (new plus literature respectively) PDLAs. 

In addition to the column densities available in the litera- 
ture (Table [8j we determine the column density of SII towards 
J2340-00 from the extant HIRES spectrum. All 3 transitions in 
the SII triplet at AAA 1250, 1253, 1259 A are well detected. We de- 
termine the column densities using the AODM and find that all 3 
lines yield values in excellent agreement. 
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QSO 


Zem lit. 


-obs 


AVlit. (kms^ 1 ) 


Zem new 


AV new (km s ) 


Include? 


1 1 .J / TV/It 


1 .990 


1.944 


4651 


1 QQ70+0 003 


4851 


N 


T9^4.0— on 


2 090 


2 054 


3516 




2825 


Y 


07222-3939 


2.18 


2.154 


2463 


7 1 R'52+O 007 


2765 


Y 


S91A 


7 7S 


7 774 


741 n 

Z4 1U 




7^QO 


IN 


Q0841+12 


2.500 


2.476 


2064 


2.4934±0.003 


1498 


Y 


B0405-331 


2.570 


2.569 


84 


2.5775±0.006 


714 


Y 


Q2343-BX415 


2.57393 


2.5720 


162 


2.5742±0.0005 


22 


Y 


Q0528-2505 


2.779 


2.812 


-2608 


2.7783±0.007 


-2664 


Y 


Q 1354- 1046 


3.007 


2.967 


3010 


3.0112±0.010 


3324 


N 


Q2O59-360 


3.09 


3.083 


514 


3.0974±0.009 


1056 


Y 


SDSS2100-0641 


3.140 


3.0924 


3469 


3.1295 


2707 


Y 


J0900+42 


3.290 


3.2458 


3107 


3.2954 


3484 


N 


J0255+00 


3.988 


3.915 


4423 


3.9936 


4759 


N 


PSSJ0957+33 


4.227 


4.178 


2826 


4.2088 


1779 


Y 



Table 7. Candidate PDLAs taken from the literature. Redshifts without error bars are taken from Prochaska, Hennawi & Herbert-Fort (2008). The final column 
indicates whether the absorber is included in our PDLA sample, which requires that the new be AV < 3000 km s . 



QSO 


Zabs 


log N(HI) 


log N(FeII) 


log N(ZnII) 


log N(SiII) 


log N(CrII) 


log N(SII) 


[Si/H] 


Reference 


J2340-00 


2.054 


20.35±0.15 


> 14.97 


12.63*0.08 


15.17*0.04 


<12.90 


14.94*0.05 


-0.69 


12,13 


Q2222-3939 


2.154 


20.85±0.10 


14.42±0.03 


<11.70 


14.55*0.05 


12.77*0.04 


14.08*0.02 


-1.81 


4,5 


Q0841 + 12 


2.476 


20.78±0.08 


14.50±0.03 


11.69*0.10 


14.99*0.03 


12.89*0.06 


14.48*0.10 


-1.30 


6 


B0405-331 


2.569 


20.60±0.10 


14.33 


<12.74 


14.74 


<13.32 




-1.37 


1 


Q2343-BX415 


2.5720 


20.98±0.05 


15.24±0.02 


12.90*0.06 


15.79*0.04 


13.58*0.04 


15.38*0.03 


-0.70 


11 


Q0528-2505 


2.812 


21.11*0.04 


15.47*0.02 


13.27*0.03 


16.01*0.03 


13.65*0.12 


15.56*0.02 


-0.61 


2,3 


Q2059-360 


3.083 


20.98±0.08 


14.52±0.07 




14.80*0.05 




14.41*0.04 


-1.69 


7 


SDSS2100-0641 


3.0924 


21.05±0.15 


15.36±0.03 


<13.14 


15.89*0.02 


13.69*0.04 


15.52*0.01 


-0.67 


8 


PSSJ0957+33 


4.178 


20.65±0.15 


14.13*0.05 




14.56*0.01 




14.39*0.06 


-1.60 


9,10 



Table 8. Column densities (in cm 2 ) for final literature PDLA sample. References: 1: Akerman et al. (2005); 2: Lu et al. (1996); 3: Centurion et al. (2003); 
4: Noterdaeme et al. (2008); 5: P. Noterdaeme private communication; 6: Dessauges-Zavadsky et al. (2007); 7 Srianand et al. (2005); 8: Herbert-Fort et al. 
(2006); 9: Prochaska et al. (2001); 10: Prochaska et al. (2003b); 1 1 : Rix et al. (2007); 12: Prochaska et al. (2007b); 13: this work. 



5.2 DLA comparison sample 

To compare the PDLAs with intervening DLAs we used the sample 
of Dessauges-Zavadsky et al. (in prep.). In order to circumvent the 
large uncertainties in AV incurred through z cra measurements, we 
impose a lower limit of 10,000 km s -1 which yields a sample of 
180 intervening DLAs with log A'(Hi) > 20.3 from the Dessauge- 
Zavadsky catalogue. We adopt the column density measurements 
compiled by Dessauges-Zavadsky et al., and convert to abundances 
using the same solar scale that was applied to the PDLAs (Asplund 
et al. 2005). In Sections [6] and [7] we compare the properties of the 
PDLAs with the DLAs in order to gain insight into whether they 
represent similar populations and what the effect of QSO proximity 
is. 



6 IONIZATION 

The proximity of PDLAs to their background QSO naturally leads 
to the question of whether they are strongly affected by the quasar's 
intense ionizing radiation. Thus far, there are conflicting indications 
in the literature for studies of individual systems. Rix et al. (2007), 
in their detailed study of a single PDLA towards Q2343-BX415, 
concluded that ionization corrections were small. Nonetheless, they 
identify absorption from highly ionized species such as NV which 
are relatively rarely detected in intervening absorbers and require 
high energy photons if produced via photoionization. One of the 



PDLAs in our sample (J1240+1455) has previously been reported 
to exhibit NV in the low resolution SDSS spectrum (Hennawi et al. 
2009) which we confirm in our echelle spectra (see Section |4"31 >. 
Fox et al. (2009) find that ~ 13% of intervening systems also ex- 
hibit NV and do not find any increase in its incidence in PDLAs 
with relative velocities 500< AV < 5000 km s~' relative to the in- 
tervening systems at A V > 5000 km s _1 . In a general survey of NV 
at all A(H 1) column densities, Fechner & Richter (2009) do find a 
proximate excess, a conflict which appears to be linked to the ve- 
locity range used (see discussion below). Finally, Lyo- emission in 
the DLA trough has been found for several PDLAs (M0ller & War- 
ren 1993; M0ller et al. 1998; Leibundgut & Robertson 1999; Elli- 
son et al. 2002; Hennawi et al. 2009), a feature which is much rarer 
amongst the intervening DLA population. In this section, we assess 
various indications of ionization in the sample of 16 PDLAs, pay- 
ing special attention to elements which might indicate ionization by 
a hard radiation source, as might be expected from the proximity to 
the QSO. 



6.1 Aluminium 

A natural probe of the level of ionization in a DLA is to compare 
column densities of different ionization states of a given element. 
The relative column densities depend on the ionization potentials, 
the shape of the ionizing spectrum and the density of the gas (or 
equivalently, the ionization parameter). The ratio of A1III/A1II is 



The nature of proximate DLAs 15 



1.5 

x 1 

E 0.5 



x 1 
E 0.5 

i.9 
x i 

3 

E 0.5 

1.8 



1 



0.5 

1.8 



X 



1 

0.5 



1.3 

x 1 
E 0.5 




CIV 1550 




: l , l\v K4d : 



— I 1 1 1 1 1 1 1 1 1 1 1 p H 

NV 1233 : 



-1000 -500 

Velocity (km/s) 



Figure 7. CIV, OVI and NV coverage in the PDLA towards Jl 131+6044 
over an extended velocity scale. The lower solid (red) line shows the er- 
ror array. OVI is detected at v ~ -750 km s _I and CIV is present from 
-1150 to approximately -500 km s -1 . There is no definitive detection of 
NV. Velocities are plotted relative to z a t> s = 2.87562. 
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Figure 12. A1III/A1II (top) and AHII/Sill (bottom) as a function of N(Hi) 
for intervening DLAs (open points, AV > 10000 km s~' ) and PDLAs (filled 
points, AV < 3000 km s -1 ). 



sometimes used as a crude estimate of ionization, although Sill is 
often substituted for A1II which saturates quickly. Vladilo et al. 
(2001) showed a broad anti-correlation of A1III/A1II with N(Hi) 
supporting the idea that this ratio traces the amount of ionized gas. 
In Figure [12] we show A1III/A1II and AHII/Sill for our PDLA and 
comparison DLA sample, where the latter ratio is corrected for the 
solar value of (Si/Al). There is no obvious difference between the 
DLA and PDLA ratios, certainly there is no indication of elevated 
amounts of A1III at a given A'(Hi) in PDLAs relative to DLAs. 
In fact, several PDLAs appear to have fairly low AHII/Sill. How- 
ever, as discussed by Howk & Sembach (1999) A1III may actually 
be a fairly poor tracer of ionized gas in intervening absorbers and 
may even originate from a different region in the absorber (Vladilo 
et al. 2001). Moreover, the behaviour of A1III/A1II is very sensi- 
tive to the shape of the ionizing spectrum. Vladilo et al. (2001) 
use Cloudy photoionization models to demonstrate that for a hard 
(QSO) ionizing spectrum the observed A1III/A1II is essentially flat 
as a function of N(H i), for a given ionization parameter (log U). 
Therefore, whilst adopting a stellar ionizing spectrum seems to nat- 
urally reproduce the trend of decreasing A1III/A1II with increasing 
,/V(H i) for a given ionization parameter, assuming a hard spectrum 
requires the log U to decrease at higher /V(Hi) to produce the same 
anti-correlation (see also Vladilo et al. 2003). It requires a very in- 
tense power law spectrum (log U > - 1) in order to achieve N(AIIII) 



> N(A1II), so the relatively low ratios of AHII/Sill observed in Fig- 
ure Q2] do not rule out proximity to a hard ionizing source. In fact, 
a hard ionizing spectrum may actually be the reason for the lower 
AHII/Sill ratios at log N(H i) <21 in the PDLAs compared to inter- 
vening DLAs. The stellar ionizing models of Vladilo et al. (2001) 
generally predict higher A1III/A1II at low A'(Hi) for soft spectra 
relative to QSO-like radiation. 



6.2 Silicon 

Next, we consider the relative abundances of Sill and SilV. Sill is 
detected in all 7 PDLAs in our sample. There are 4 SilV detec- 
tions, two lower limits (from saturated lines) and one upper limit 
from a non-detection. For the two PDLAs with log ,/V(H i) < 20.4 
(Q0151+048 and J0142+0023) we find that log[N(SiIV)/N(SiII)] = 
-0.26 and -0.42 respectively. Such high fractions of SilV are rel- 
atively rare in intervening DLAs. The other PDLAs where SilV 
has a detection or upper limit have SilV/Sill fractions less than 
10%. In cases where intervening DLAs exhibit SilV, it is usu- 
ally observed in a different velocity structure from Sill (Wolfe & 
Prochaska 2000), which can be obtained if the ionizing spectrum 
is relatively soft (Howk & Sembach 1999). This is indeed the case 
for the PDLAs towards J2321 + 1421 and J0142+0023 where the 
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bulk of the SilV is offset from the Sill. However, for 3/6 of the PD- 
LAs in our sample where we detect SilV, it traces the structure of 
the Sill (Q0151+048, J1131+6044, J1604+3951). This is unusual; 
one explanation could be the presence of a hard spectrum (and/or 
high ionzation parameter). For the final PDLA where SilV is de- 
tected (J1240+1455) the SilV is so strongly saturated we have no 
information on its velocity structure. Interestingly, the coincident 
(or not) velocity structure of SilV with Sill is independent of the 
AV of the PDLA. For example, J2321 + 1421 has the smallest ve- 
locity separation (-1616 km s -1 ) with little or no SilV coinciding 
with Sill, yet the largest AV PDLA (Jl 13 1+6044, 2424 km s" 1 ) has 
very similar velocity structure in the two ions. 

6.3 Alpha elements 

S, Si, Ar and O are all alpha capture elements that are produced 
predominantly in massive stars. Their shared nucleosynthetic ori- 
gin leads to approximately solar abundance ratios of these elements 
in Galactic stars (e.g. Chen et al. 2002). Departures from solar ra- 
tios in DLAs can occur when one or both of the alpha elements 
in question is depleted from the gas phase onto dust, or if there is 
significant ionization. 01 is an ideal 'anchor' to study alpha ele- 
ments since oxygen is both relatively undepleted (Savage & Sem- 
bach 1991) and 01 requires a negligible ionization correction. Un- 
fortunately, the most accessible 01 transition (A = 1302 A) is usu- 
ally strongly saturated. We have only one PDLA in our sample with 
a well-constrained column density of OI; the other absorbers yield 
only lower limits. 

We have detections or meanginful limits for S and Si for all the 
PDLAs in our sample. Like O, S is largely undepleted in the Galac- 
tic ISM, but Si is mildly refractory. Super-solar ratios of [S/Si] may 
therefore be observed if dust is present. In Figure[T3lwe plot [S/Si] 
versus A^H i) for the sample of 7 new PDLAs plus systems taken 
from the literature (see Table[8]l. PDLAs are colour-coded by their 
relative velocity and intervening DLAs with AV > 10,000 km s -1 
are shown for comparison as open circles. The intervening DLAs 
cluster around [S/Si]=0 with most points within ±0.2 dex of the so- 
lar ratio. The effects of dust depletion appear to be relatively mild 
(i.e. few DLAs show highly super-solar [S/Si]). The PDLAs also 
cluster around the solar ratio, but only when log A^H i) > 20.75. At 
least 50% (accounting for limits) of the lower A'(Hi) PDLAs have 
[S/Si] that is several tenths of a dex below that of the intervening 
systems. Since neither nucleosynthetic nor depletion are expected 
to lower the [S/Si], the deviation from the solar ratio may be due 
to ionization effects. This hypothesis is supported by the tendency 
of low [S/Si] ratios to be observed only at lower iV(H i) PDLAs. 
Photoionization modelling indeed indicates that S and Si require 
corrections of opposite sign in the presence of a power law spec- 
trum (Rix et al. 2007). Interestingly, the trend of low [S/Si] at low 
A'(Hi) does not seem to depend solely on AV. Even at AV > 1500 
km s -1 , 2/3 PDLAs with log A'(Hi) < 20.5 show very sub-solar 
[S/Si] ratios. Conversely, PDLAs with very small velocity separa- 
tions exhibit solar [S/Si] when log N(Hi) > 21, indicating that they 
are well-shielded despite a small velocity separation. However, we 
caution that many of the relative velocities are constrained only to 
within ± 500 km s and that AV is not necessarily an indication of 
physical proximity. 

A natural concern in the interpretation of Figure [13] is that it 
relies on several limits in PDLAs resulting from the non-detection 
of SII. As discussed in Section [4] the calculated limit is somewhat 
sensitive to the assumed width of the line. In order to check our 
calculated SII detection limits we produced Voigt profiles of SII 
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Figure 13. [S/Si] as a function of log N(H i) for intervening DLAs (open 
points, AV > 10000 km s~') and PDLAs (filled coloured points) from 
the new data presented here and taken from the literature (see Table [8). 
PDLAs are colour coded by their relative velocities (see Table ff): cyan 
for AV < 500 km s -1 , green for 500 < AV < 1500 km s _1 and red for 
AV > 1500 km s _I . The low values of [S/Si] can not be easily explained 
by dust depletion or nucleosynthesis and may be caused by ionization by a 
hard radiation source. 

A 1259 for the limiting column density, adopting the b- value and 
redshift of the strongest component in the fit to one of the low ions 
(either Sill or Fell). A visual inspection of the calculated limits 
overlaid on the data indicate that they are reasonable estimates. 

Argon is relatively rarely measured in DLAs, mostly due 
frequent blending of the rest-frame far UV lines of Arl at A = 
1048, 1066 A. However, the abundance of Arl is very sensitive 
to ionization, because its ratio of photoionization to recombina- 
tion rates are typically one order of magnitude higher than for HI 
(Sofia & Jenkins 1998). Interestingly for our study, Vladilo et al. 
(2003) have used Cloudy modelling to demonstrate that the frac- 
tion of Arl is sensitive to the adopted radiation field. Low ratios 
of Ar relative to other elements are expected when the radiation 
field is hard, whereas the models predict solar ratios when the ion- 
izing spectrum is soft. Two of our limits on [Ar/Si] are not deep 
enough to be very meangingful: [Ar/Si] < +0.42,-0.01 towards 
J0140-0839 and J2321 + 1421 respectively. The two other limits 
are very sub-solar: [Ar/Si] < -0.86,-0.47 towards Jl 131-6044 
and J0142+0023 respectively. We have one new argon detection 
(towards J1604+3951), although we can only derive its ratio with 
silicon for the principal component: [Ar/Si] =-0.40. Sub-solar ra- 
tios of [Ar/Si] are therefore apparently common in the PDLAs and 
further support significant ionization by a hard radiation field. 

6.4 Nitrogen 

The most highly ionized species available for study in our sample 
is NV (in all cases OVI is blended). Fox et al. (2009) find a 13^% 
NV detection rate in DLAs compared with 13*i 8 % in PDLAs out to 
5000 km s from the QSO. This result may seem surprising given 
the ionization radiation from the nearby QSO. Indeed, Fechner & 
Richter (2009) find a high NV incidence rate amongst AV < 5000 
km s" 1 absorbers with 13 < log N(Hi) < 17.0. Fox et al. (2009) 
suggest that the consistency in NV detection rates between DLAs 
and PDLAs may be due to the fairly large velocity interval over 
which they include a DLA in their proximate sample. In our sam- 
ple, 2/7 (J1240+1455 and J1604+3951) of the PDLAs exhibit NV 
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at a similar velocity to the low ions - tentative evidence that NV 
is more common in PDLAs, although better statistics are required 
to confirm this. Notably, the PDLA towards J1240+1455 exhibits 
one of the largest N(NV) yet reported in the literature (Fox et al. 
2009). The two NV detections in our sample occur in PDLAs with 
small AV, high N(H i) and relatively high metallicity (Z ~ 1 / 10Z G ). 
Rix et al. (2007) also detected NV in the iV(Hi)=20.98 PDLA to- 
wards Q2343-BX415, which also has a fairly low velocity separa- 
tion (AV = 22 km s -1 ) and high metallicity (Z ~ 1 /5Z Q ). However, 
J2321 + 1421 which has the largest negative velocity in our sample 

( 1600 km s -1 ) does not exhibit NV. The detection rate of NV 

would therefore increase if we limited our statistics to lower AV, 
supporting the explanation of Fox et al. for the lack of NV excess 
in their proximate sample. The higher incidence of NV out to 5000 
km s" 1 seen by Fechner & Richter (2009) may be associated with 
the lower jV(Hi) column densities of their sample. 




7 ABUNDANCES 

7.1 Metallicity 

Figure [14] shows Fe and Si abundances as a function of redshift for 
DLAs and PDLAs whose redshifts are Zabs > L8. Our sample in- 
cludes PDLAs that are amongst both the most metal-poor and the 
most metal-rich for their redshift. The PDLA towards JO 140-0839 
has the lowest metallicity ever detected in a high yV(H i) absorber 
(see Section |4~T| for further discussion), although there may be an 
upward correction of up to a few tenths of a dex to be made for 
ionization. The results in Figure [14] ostensibly support the conclu- 
sion of Rix et al. (2007) that PDLAs exhibit metallicities that span 
the distribution of intervening DLAs. However, in Section [6] it is 
argued that the low A(Hi) systems are likely to have significant 
ionization corrections. We therefore re-assess the PDLA metallici- 
ties as a function of A(H i). 

PDLAs with log N(H i) > 2 1 probably have negligible ioniza- 
tion corrections and Figure Q3] shows that these PDLAs have quite 
high [Si, S/H]. The Zn abundances are also distributed towards the 
upper end of the intervening DLA distribution. [Fe/H] shows more 
scatter and its sensitivity to dust depletion makes it harder to in- 
terpret. At lower A'(Hi), all but one of the PDLAs have Si and S 
abundances that are below the median. This intriguing observation 
could be interpreted as a difference in the provenance of the low 
and high Af(Hi) PDLAs. Alternatively, (and perhaps more likely) 
it could indicate that the PDLAs (or at least a subset thereof) have 
intrinsically higher metallicities, but whose abundance determina- 
tions are affected by ionization at low A(Hi). 

7.2 Dust depletion indicators 

In Figure[T6lwe plot the [Zn/Fe] ratio (as an indicator of depletion) 
as a function of [Zn/H] for the PDLA and DLA samples. The me- 
dian [Zn/Fe] for our literature sample (when both Zn and Fe are 
measured) is [Zn/Fe] — 1-0.4. J1240+1455 appears to have partic- 
ularly high depletion, although at Z ~ l/10Z o it is also one of the 
more metal-rich absorbers. The PDLA towards J1240+1455 has a 
high A'(Hi) so the high [Zn/Fe] ratio is unlikely to be caused by 
ionization. J1604+3951 shows clear signs of different levels of de- 
pletion between two velocity components with high [Zn/Fe] in one 
component but not the other (see Section |4~6t . Overall, there is no 
obvious difference between the dust depletion in the PDLAs and 
the intervening systems. 
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Figure 14. Si and Fe abundances as a function of absorption redshift for 
intervening DLAs (open points, AV > 10000 km s _1 ) and PDLAs (filled 
coloured points). PDLAs are colour coded by their relative velocities (see 
Table[8j: cyan for AV < 500 km s~\ green for 500 < AV < 1500 km s" 1 
and red for AV > 1500 km s -1 . There is a large scatter in both Si and Fe 
abundances in PDLAs at a given redshift (see Section[TT}- 



7.3 Alpha elements 

[Si/Fe] is used as an indicator of a enhancement, since both ele- 
ments are measured in the majority of our sample. [Si/Fe] is usually 
plotted as a fuction of metallicity, although the results presented 
so far in this paper indicate that trends may also be present with 
A(Hi) due to ionization. In Figure [T7] we therefore plot [Si/Fe] as 
a function of both Af(Hi) and [Si/H]. Figure [171 shows that at low 
A(H i) the PDLAs have fairly typical values of [Si/Fe], with the ex- 
ception of the very metal-poor PDLA towards J0140-0839 which 
has [Si/Fe] > 0.72. This PDLA is even more of an outlier when its 
[Si/Fe] is considered relative to its metallicity. Although the high 
[Si/Fe] ratio could be due in part to ionization corrections and/or 
dust, it also raises the intriguing possibility of extreme alpha ele- 
ment enhancement in a chemically young object. SDSS offers the 
opportunity to search out other rare, low metallicity DLAs in order 
to investigate their nucleosynthetic histories. The [O/Fe] ratio in 
this absorber is also high, [O/Fe]^ 0.75 and OI is much less affected 
by ionization and dust than Sill. Jl 13 1+6044 also has a relatively 
high value of [Si/Fe]=0.67, although Figure [T7] shows that such 
values are not unknown in the intervening population. J1240+1455 
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Figure 15. Si, Fe, Zn and S abundances as a function of log A'(Hi) for intervening DLAs (open points, AV > 10000 km s ) and PDLAs (filled coloured 
points). Si and Fe are either mildly or severely depleted by dust, whereas Zn and S are largely undepleted from the gas phase and may therefore represent more 
reliable indicators of metallicity. PDLAs are colour coded by their relative velocities (see Table[8j: cyan for AV < 500 km s , green for 500 < AV < 1500 
km s and red for AV > 1500 km s . Silicon and sulphur (and to a lesser extent, zinc) exhibit high abundances when log A'(Hi) > 21 (see Section PTTV 



has an extremely high [Si/Fe] that is unlikely to be due to ionization 
effects, given its high A(H i). However, as discussed in the previous 
subsection, the [Zn/Fe] indicates a large depletion fraction which 
could lead to an over-estimate of [Si/Fe] (Fe is much more refrac- 
tory than Si). Using the undepleted ratio of [S/Zn]=0. 11, the high 
[Si/Fe] of J1240+1455 appears to be largely due to dust. In sum- 
mary, there are no clear, systematic differences in the [Si/Fe] ratios 
of PDLAs compared with the intervening systems, although some 
PDLAs do appear to have relatively high values. 



8 CONSTRAINING THE DISTANCES TO PDLAS FROM 
FINE STRUCTURE LINES 

We have emphasized in this paper that the velocity offsets (AV) be- 
tween the QSO and PDLAs are unlikely to be convertible into dis- 
tances using Hubble's law. To attempt to constrain the true physical 
distances we can model the interplay between the QSO's radiation 
and the PDLA's interstellar gas. 



8.1 Distances from CII* 

In this subsection we describe the method for using C II* to esti- 
mate the PDLA-QSO distance, as well as describing the assump- 



tions of the model, limitations and caveats. The practical applica- 
tion of the method is given in Section [8.1.1l 

Wolfe, Prochaska & Gawiser (2003) have used CII* to esti- 
mate the radiation intensity, and hence the star formation rate in 
DLAs. The far UV radiation emitted by massive stars contributes 
to heating through the grain photo-electric effect, which in turn 
heats the gas. The heating rate, T is a function of the mean radi- 
ation intensity (J), the heating efficiency (e) and the dust to gas 
ratio relative to the Galactic value (k): T cc Jck. The value of e is 
known from Galactic studies (e.g. Bakes & Tielens 1994), and k 
can be inferred from ratios of refractory to undepleted elements for 
individual DLAs, under an assumed model for dust depletion pat- 
terns and intrinsic abundance ratios. In a plane parallel layer J is 
proportional to if/* , the star formation rate per unit area (which we 
henceforth refer to as simply 'star formation rate' for brevity). The 
star formation rate can therefore be determined once T is known. 
This is achieved by assuming that the medium is in local thermal 
equilibrium and that heating rate can thus be equated to the cooling 
rate. This latter quantity is inferred from the C II* absorption line 
which measures the population of the excited 2 P3/2 fine-structure 
state that spontaneously decays to the ground-state by emitting a 
158//m photon. 

For PDLAs, a second source of far UV photons potentially 
contributes to the heating: radiation directly from the proximate 
QSO. Rix et al. (2007) measured CII* in one PDLA and by assum- 
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[Zn/H] 

Figure 16. [Zn/Fe] as a function of [Zn/H] for intervening DLAs (open 
points, AV > 10000 km s -1 , limits excluded for clarity) and PDLAs (filled 
coloured points). PDLAs are colour coded by their relative velocities (see 
Tabled): cyan for AV < 500 km s -1 , green for 500 < AV < 1500 km s -1 
and red for AV > 1500 km s . The PDLAs show no distinct dust depletion 
trends compared with the intervening absorbers (see Section l77H . 



ing that all the radiation originated from the QSO placed a lower 
limit on the distance of the PDLA from the QSO. We follow a sim- 
ilar procedure here, but calculate the distance as a function of i/r* 
to explore the range of likely values of the PDLA-QSO separation. 
We assume that the mean radiation intensity, J, inferred from CII* 
has contributions from star formation and the QSO: 



TOT - J&F + JqSO 



(5) 



The radiation contribution from the QSO depends on the QSO lu- 
minosity (L v ) and separation (d). J is measured in ergs ctrT 2 s -1 sr~' 
Hz -1 so that 

Ly 

And 1 ' 



4nJ QSOv 



(6) 



The CII* analysis of Wolfe et al. (2003) yields J for a rest wave- 
length of 1500 A. 

Combining equations[5]and[6]and solving for d, we obtain 



d-±J — 

4n V 7tot,v - /sf.v 



(7) 



The contribution to I from star formation depends on the geometry 
of the plane parallel layer (parameterised by the ratio of its radius 
and height, R/h), metallicity and its dust-to-gas ratio, k v . For k v h <K 
k v R <K 1 equations 18 and 19 of Wolfe et al. (2003) can be reduced 
and re- arranged to give 



Jsf.v = 8.4 x lO-'V* 



1 + ln(R/h) 
Sn 



(8) 



where i/r* is measured in units of Mq yr~' kpc~ 2 . A fixed aspect 
ratio of R/h=20 is used for all calculations. Following Wolfe et al. 
(2003), the dust-to-gas ratio is determined from Si/Fe and metallic- 
ity from [Si/H] with an assumed intrinsic [Si/Fe] =0.2 and SMC de- 
pletion patterns. Solutions can be calculated for two cases: gas that 
is dominated by the cold and warm neutral media (CNM, WNM) 
respectively. However, using CII* to constrain distances for the 
WNM solutions has two problems. First, in the WNM, the cool- 
ing is actually dominated not by the [CII] 158 yum transition, but by 
Lya and the recombination of electrons onto grains. Second, at the 
low densities implied by the WNM solutions, the dominant heating 
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Figure 17. [Si/Fe] as a function of [Si/H] (bottom panel) and log N(H i) (top 
panel) for intervening DLAs (open points, AV > 10000km s~') and PDLAs 
(filled coloured points). PDLAs are colour coded by their relative velocities 
(see Table[8j: cyan for AV < 500 km s~\ green for 500 < AV < 1500 
km s _1 and red for AV > 1500 km s~' . The PDLAs show no distinct [a/Fe] 
trends compared with the intervening absorbers (see Section[73}. 



source is no longer the grain photo-electric effect, but cosmic ray 
heating and, to a lesser extent, heating by X-rays. In the calculations 
of Wolfe et al. (2003), the cosmic ray production rate is assumed 
to scale with the star formation rate. However, we are additionally 
concerned with energy sources associated with the QSO and it is 
unclear whether AGN produce cosmic rays and at what rate. We 
therefore only consider the CNM solutions, but caution that these 
are probably inappropriate for the PDLAs with the lowest cooling 
rates. A full treatment of the WNM solutions for AGN heating is 
beyond the scope of this paper. 

The QSO luminosity at Aq = 1500 A is determined from the 
SDSS spectrum, except in the case of Q0151+048. For this target, 
we take B AB = 17.83 from Fynbo et al. (2000) which converts to a 
flux of 2.7xl0~ 27 ergs s~' Hz~' (no corrections are made for hya 
emission or Lyo- forest absorption, since they lie outside the B filter 
for z = 1.9). The luminosity and flux (in units of ergs s~' Hz~') at 
a given frequency are related through the equation 



F v = L, 



(1+Z)V ' 



l+z 
~An~d\ 



(9) 
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Figure 18. Distances between the PDLA and Q015 1+048 calculated from 
the CII* analysis. The curve shows the distance corresponding to an in- 
creasing local contribution of radiation from star formation rate, log if/*. 
The vertical dotted line indicates log i//* = -2.4 M Q yr . kpc 




G/G 

Figure 19. Ratio of fine structure to ground state column densities from 
UV pumping determined from the popratio software. The arrows indicate 
the inferred upper limit of the incident radiation field in units of the Galactic 
value based on upper limits of N(SiII*). 



where d L is the luminosity distance in our adopted cosmology. In 
Table|9]we list the calculated values of Lisoo- 

In previous sections, we have argued that ionization may lead 
to errors in in the column densities that we measure, particularly 
at low 7V(Hi). The metallicity (nominally [Si/H]) and dust to gas 
ratio (calculated from [Si/Fe]) are required to calculate J T ot- How- 
ever, Wolfe et al. (2003) have shown that J T0T is actually relatively 
insensitive to the metallicity used for values Z < l/10Z o . Even 
at the highest metallicities, incorrect values of [Si/H] are likely to 
affect J tot by less than 20%. Uncertainties in the dust-to-gas ra- 
tio are potentially more problematic, since the grain photo-electric 
effect can be a dominant source of ISM heating. However, the con- 
tribution of this process to the heating budget depends on the ISM 
density (e.g. Figure 3 in Wolfe et al. 2003), becoming more domi- 
nant at high densities. For CNM solutions, lower dust-to-gas ratios 
(i.e. fewer grain targets for heating) requires a higher incident in- 
tensity to account for the observed CII*Ll However, in Section 1731 
we showed that the range of [Si/Fe] ratios of most of the PDLAs are 
in good agreement with intervening systems. Assuming that there 
is an intrinsic 'floor' to [Si/Fe] whose value is approximately 0.2 
(e.g. Prochaska & Wolfe 2002) and given that ionization effects 
tend to lead to over-estimates of [Si/Fe] (e.g. Dessauges-Zavadsky 
et al. 2003; Howk & Sembach 1999), the measured Si/Fe ratios do 
not indicate significant corrections in the majority of cases. 

8.1.1 Application to the data 



CII* is detected in Q0151+048 and J1604+3951 and we 
have upper limits for J0140-0839, Jl 131+6044 and J2321 + 1421 
which will yield lower distance limits for a given i/r*. However, 
the CII* towards J1604+3951 is partly blended with CII A 1334, 
and the unblended components are saturated, so we do not consider 
this PDLA further. Q015 1+048 may have ,/qso contributions from 
Q0151+048B, introducing a second unknown distance. Fynbo et al. 

7 As discussed in the text, at low densities cosmic ray heating dominates 
over the grain photo-electric effect. WNM solutions, although not pursued 
here, are therefore very insensitive to the input value of [Si/Fe]. 



(1999) have argued that the most likely orientation of the system 
is that Q0151+048B is in front of Q015 1+048 and the DLA and 
causes Lya- emission on the near face of the absorber. Modelling 
this system completely is complicated not just by uncertainty in the 
geometry of the QSO pair and absorber, but also by the anisotropy 
of quasar radiation and hence the flux 'seen' by the DLA com- 
ing from Q0151+048B. However, assuming that all the flux comes 
from the background QSO still yields a lower limit to the distance 
between Q0151+048 and the PDLA (just as we determine a lower 
limit under the assumption that internal star formation does not 
contribute). 

In Table [9] we list the limiting value of the PDLA-QSO dis- 
tance for a star formation rate if/* = and also the fiducial case of 
a Galactic star formation rate of log if/* = -2.4 M Q yr~' kpc ~ 2 for 
the CNM solutions. For PDLAs with upper limits to N(CII*) the 
CNM solutions are likely to be inappropriate and they are given 
in the Table only for completeness. More sensitive CII* limits 
could definitively rule out CNM solutions entirely as is the case for 
J1131+6044 which only has a WNM solution. The CII* analysis 
therefore only yields useful a distance constraint for Q0151+048 
where the lower limit is 160 kpc (no star formation case), rising to 
just over 400 kpc for a Galactic star formation rate, a separation 
approximately 5 times larger than the only other PDLA in the lit- 
erature has been analysed in this way (Q2343-BX415 by Rix et 
al. 2007). In Figure [T8] we show the complete range PDLA-QSO 
distances as a function of star formation rate for Q0151+048. The 
Figure demonstrates that once star formation dominates over the 
QSO's radiation, the inferred distance to the PDLA rises dramati- 
cally, essentially providing an upper limit to the likely star forma- 
tion rate in the galaxy, as well limits of the distance to the QSO. 



8.2 Distance constraints from Sill* limits 

There are a number of uncertainties and assumptions that underpin 
the distance calculations from CII* which affect the inferred radi- 
ation field to varying extents. These include knowledge of the gas 
phase (cold or warm), metallicity, dust depletion, assumed geom- 
etry and extinction law. An alternative assessment of the distance 
between the absorber and the radiation source can be made using 
the Sill fine structure lines and the ratio of N(SiII*) to the column 
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QSO 




AV 


log N(CII*) 


log N(SiII*) 


Ll500 


dcNM (kpc) 


dcNM (kpc) 








(kms- 1 ) 


(cm- 2 ) 


(cm- 2 ) 


(ergs s Hz -1 ) 


^* =0 


log <K = -2.4 


(kpc) 


J0140-0839 


3.6960 


1250 


<12.41 


<11.48 


1.08xl0 32 


>2099.0 




>14.8 


J0142+0023 


3.34765 


1772 




<11.17 


4.76xl0 31 






>29.8 


Q015 1+048 


1.9342 


-1199 


13.0 


< 11.81 


2.44xl0 31 


160.5 


406.0 


>8.6 


J1131+6044 


2.8754 


2424 


<12.51 


<1 1.57 


5.85X10 31 






>30.8 


Q2321+142I 


2.5731 


-1616 


<12.55 


<11.78 


2.27X10 31 


>510.1 




>14.3 



Table 9. PDLA-QSO distance calculations. Solutions are given for the CII* analysis under the assumption of a CNM and for Sill* where suitable detec- 
tions/limits are possible (see text for details). 



of the ground state ions. In order to determine the relation between 
N(SiII*)/N(SiII) and the incident radiation field, we use the pub- 
licly available code popratio (Silva & Viegas 2002). We assume 
that Sill* is populated purely by UV-pumping (fluorescence) and 
ignore contributions from direct (IR) excitation and collisions. The 
results from popratio are shown by the curve in Figure[l9]as a func- 
tion of the radiation field in Galactic units (Go). Upper limits for 
N(SiII*) are derived from Sill*/! 1264 which in turn yield upper 
limits on G/Gq, which is related to J 

J = G x 1 x 10-' 9 ergs cnT 2 s~' sf" 1 Hz" 1 . (10) 

The luminosities in Table[9]are then combined with J in equa- 
tion[6]to give a lower limit on the distance. These values are lower 
limits not just because of the non-detection of Sill*, but also be- 
cause of our assumption that only UV-pumping contributes to the 
population of the fine structure level. The lower limits derived are 
typically 15 to 30 kpc, see Table [9] Although these are consider- 
ably less stringent than the CII* distances, they are fairly robust 
and depend on very few assumptions and input parameters. The 
only caveat to this analysis is that at least two of the Sill transitions 
from the ground-state must be optically thin (Sarazin, Rybicki & 
Flannery 1979), which holds for all but one of our PDLA sample 
(which is excluded from this analysis). 

In summary, we have used two methods to constrain QSO- 
PDLA distances. The CII* method (under the CNM assumption) 
gives a lower limit of 160 kpc for Q0151+048 (assuming no in- 
ternal star formation), but is rather dependent on model assump- 
tions. The Sill* model is more robust but gives less stringent lim- 
its, typically > 15-30 kpc. Adding star formation to either method 
increases the inferred distance of the PDLA from the QSO. The 
distances imply that the absorbers are external to the QSO. 



9 DISCUSSION 

9.1 Sulphur and argon as indicators of QSO proximity 

We have suggested that sub-solar ratios of [S/Si] at low values of 
A'(Hi) may be caused by the significant under-estimate of N(S) 
from SII and, to a lesser extent, an over-estimate of N(Si) from 
Sill in the presence of a hard ionizing spectrum. This result could 
be predicted from the photoionization models of Rix et al. (2007), 
who showed that the observed under-estimate of sulphur is strongly 
dependent on the ionization parameter, U. Most models of interven- 
ing DLAs and sub-DLAs have concluded that the value of log U is 
typically < -3 (e.g Howk & Sembach 1999; Dessauges-Zavadsky 
et al. 2003). However, for absorbers close to a QSO, the situation 
can be quite different. For example, Prochaska & Hennawi (2009) 
find that for a hydrogen volume density of 0.1 atoms crrT 3 the ion- 
ization parameter of one of their transverse sub-DLAs at a distance 



of ~ 100 h- l Q kpc from the QSO is log U = -1.5. The QSOs in our 
sample are typically a factor of 10 brighter than the case studied 
by Prochaska & Hennawi (2009), so ionization parameters in the 
range -2 < log U < are quite feasible. 

Sub-solar [S/Si] ratios may therefore be used as a signpost of a 
nearby hard radiation source. Two of the intervening DLAs (FBQS 
2334-0908 and PSS 0133+0400) in our literature sample also have 
[S/Si]< -0.5. One possible explanation could be that although they 
are intervening (AV » 10,000 km s~') systems, there is a nearby 
foreground QSO at a similar redshift. FBQS 2334-0908 is covered 
by the SDSS footprint. Although none of the objects near to FBQS 
2334-0908 were targeted by the SDSS for spectroscopy, there are 
two point sources at separations of 27 and 105 arcseconds (208 and 
809 hzl kpc at z = 3 respectively) whose colours are consistent 
with expectations of QSOs. PSS 0133+0400 is not covered by the 
SDSS, but the APM catalogue shows 2 point sources with similar 
colours to the QSO at separations of 57 and 64 arcseconds (409 
and 459 kpc at z = 3.7 respectively). It would be interesting 
to obtain spectra of these sources to determine whether or not they 
are indeed QSOs at the same redshift as the intervening DLAs. In- 
deed, some of the intervening NV absorbers identified by Fox et 
al. (2009) have nearby QSOs at the same redshifts (G. Worseck, 
private communication). 

The low ratios of [Ar/Si] (where measured) in our PDLAs are 
also indicative of a hard ionizing spectrum (Vladilo et al. 2003). 
Using a sample of 10 Arl measurements, Vladilo et al. (2003) have 
also found low ([Ar/Si]< -0.5) at z abs < 3. At higher redshifts, 
the values are close to, or at, the solar level. Vladilo et al. (2003) 
have argued that this may be evidence for evolution in the ioniz- 
ing background, which is shifting from a softer to harder shape as 
the redshift decreases. However, 3/7 of the z a b s < 3 absorbers in 
the Vladilo et al. (2003) sample are PDLAs. A further two have 
velocities within 10,000 km s" 1 , leaving only two which might be 
considered truly intervening absorbers. At z a bs > 3 all 3 absorbers 
have AV > 15,000 km s _1 . Given the evidence presented in this 
paper that the effects of a hard ionizing spectrum are seen out to at 
least 3000 km s" 1 (and previous results that indicate narrow associ- 
ated absorbers may contribute out to 10,000 km s~', e.g. Wild et al. 
2008; Tytler et al. 2009), the apparent redshift evolution may actu- 
ally arise from the inclusion of proximate systems. A larger sample 
of large AV DLAs with argon measurements is necessary to explore 
more fully the redshift evolution of [Ar/Si]. 



9.2 NV as a tracer of high ionization gas 

High ionization species such as NV and OVI might be expected 
to be more common in absorbers close to the QSO. Although OVI 
is blended in almost all of our sightlines, we have detections and 
upper limits of NV for all of our PDLAs. It is intriguing that the 
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two detections of NV at the same velocity as the singly ionized 
metal lines are towards the two highest Af(Hi) absorbers (towards 
J1240+1455 and J1640+3951); both have log N(Hi)> 21. How- 
ever, both absorbers also have normal [S/Si] ratios, whereas we 
might expect departures from the solar ratio if the gas is partly ion- 
ized by a hard spectrum (as is seen for low A'(Hi) PDLAs). The 
presence of NV in an ISM that is dominated by neutral gas indi- 
cates that the NV might be formed through an internal, localised 
source of ionizing radiation. In that case, the radiation from the 
QSO may not be responsible for the formation of the NV and the 
proximate nature of the DLA is mis-leading. Fox et al. (2009) found 
that the NV detection rate is higher for higher metallicity (inter- 
vening) DLAs. J1240+1455 and J1640+3951 have (undepleted) S 
abundances that are more than 1 dex higher than many of the other 
PDLAs in our sample and 3 times higher than intervening DLAs 
at the same Af(Hi) and redshift. Fox et al. (2009) have suggested 
that active star formation could be responsible for highly ionized 
gas which is consistent with the high metallicities, large velocity 
spreads and (for J1604+3951) the very high CII* column density. 
Indeed, Lehner et al. (2008) found that a pure QSO (hard) spectrum 
always under-produced the observed amount of NV in their study 
of a DLA with multi-phase gas; adding a soft stellar component 
significantly increased the predicted NV column density. Our own 
Cloudy models indicate that even in the presence of a hard radi- 
ation field that is sufficient to produce sub-solar [S/Si] the NV/NI 
fraction is still only on the order of 10%, which is consistent with 
the NV column density limits in our lower N(H i) systems. 

Both J1240+1455 and J1640+3951 show additional (and in 
the former case, much stronger) NV offset by several hundred 
km s~' to the blue. A third PDLA (towards Q0151+048) shows 
highly offset NV (by -825 km s~') but no NV at zero velocity; in 
this case, the N(Hi) is much lower, only log N(Hi) = 20.34. It is 
notable that every PDLA in our sample of 7 echelle spectra (with 
the exception of J0140-0839, which has very weak metal lines in 
general) has high ionization gas that is offset to negative velocities. 
In many cases, there is high ionization gas (such as CIV) at the 
redshift of the low ions but no low ions are seen to accompany the 
negative velocity high ion components. In some cases the velocity 
offset is small (e.g. CIV at ~ -100 km s" 1 in J0142+0023). In other 
cases, the velocity offsets exceed -500 km s _I (e.g. Q015 1+048). 
We summarise these offset components in Table [Tol (see Sections 
14.21 to |4~7l for more detailed descriptions) listing both the AV of 
the PDLA from the background QSO and the velocity offset of the 
high ions from the main low ion components. Rix et al. (2007) also 
found a highly ionized component (observed in NV, CIV and SilV) 
offset by ~ +500 km s _I from the PDLA towards Q2343-BX415. 
Weidinger et al. (2005) find NV in 2 absorbers offset from a proxi- 
mate Lyman limit system by a few thousand km s~' . In these more 
extreme cases (e.g. Q0151+048) the absorption is likely to be un- 
connected with the proximate HI system and more akin to the as- 
sociated high ionization associated systems studied by D'Odorico 
et al. (2004) and Fechner & Richter (2009). 

For J 1 1 3 1+6044 and Q0 1 5 1 +048 where the high ion velocity 
offset is very high, it is likely that the absorption is arising in gas 
outside the PDLA. As discussed in Section 1331 Q0151+048 has a 
fainter companion separated by 27.5 /i^J kpc and — 1-120 km s . 
Although there is no confirmed companion for Jl 131+6044 in the 
NASA Extragalactic Database (NED), the SDSS image shows a 
bright (g = 17.59) point-source 5.66 arcseconds to the north-west. 
The object has an almost identical g—r colour to Jl 13 1+6044 (0.41 
and 0.38 respectively), so that this may be a binary QSO. The puta- 
tive companion may contribute to the significant ionization seen in 



the Jl 131+6044 PDLA (e.g. high SilV fraction), despite the large 
positive velocity offset from the QSO (2424 km s , but see below 
for further discussion on the interpretation of velocities). The other 
5 QSOs in our sample have no obvious bright companion within 20 
arcseconds on the SDSS image. If the companion to Jl 131+6044 
is confirmed, this would mean that 2/7 of the QSOs with PDLAs 
in our sample have a close companion. Q2343-BX415 also has a 
companion, and its PDLA also has 2 distinct velocity components 
in the highly ionized gas, one of which is coincident with the bulk 
of the low ions, and the other offset by ~ +500 km s _I . 



9.3 The origin and nature of the PDLAs 

The question that underpins the research of all associated absorbers 
is whether they are intrinsic to the QSO (host or outflow) or simply 
nearby in velocity space. M0ller et al. (1998) examined a num- 
ber of hypotheses for the nature of the PDLAs and favour a model 
in which the PDLAs are similar in nature to the intervening ab- 
sorbers, but possibly located in a preferential environment, such as 
in the same overdensity as the QSO. Rix et al. (2007) suggest that 
the PDLA towards Q2343-BX415 (included in our literature sam- 
ple) may be associated with outflowing material from the quasar 
host. Indirect clues to the provenance of the PDLAs may be gar- 
nered from their statistical properties. An excess of PDLAs relative 
to intervening systems has been confirmed by three studies (Elli- 
son et al. 2002; Russell et al. 2006; Prochaska et al 2008b). Al- 
though Ellison et al. (2002) found an additional excess of PDLAs 
in their radio-selected sample relative to optically selected QSOs, 
Russell et al. (2006) found an equal PDLA enhancement towards 
radio-loud and radio-quiet QSOs. The discrepancy may be due to 
the limited statistics of the Ellison et al. study, but a second possi- 
bility is the nature of radio sources. The CORALS sample used by 
Ellison et al. (2002) is comprised entirely of (rarer) flat-spectrum 
quasars that have compact morphologies. The sources detected at 
20-cm in the Russell et al. (2006) sample will have a range of spec- 
tral indices and orientations. It would be interesting to investigate 
the dependence of PDLA incidence as a function of radio spec- 
tral index, as has been done extensively for CIV absorbers. There 
is certainly evidence that the quasar's radiation affects the ability 
of a PDLA to survive. Hennawi & Prochaska (2007) have shown 
that the incidence of transverse DLAs in projected QSO pairs over- 
predicts the incidence of proximate absorbers by a factor of 4 - 
20. These authors suggested that line-of-sight PDLAs are preferen- 
tially photo-evaporated by the QSO's beamed radiation. A similar 
overabundance of transverse Mgll absorbers relative to the line of 
sight has been seen by Bowen et al. (2006). Moreover, Prochaska 
et al. (2008b) show that despite the higher incidence of PDLAs rel- 
ative to the intervening population, a clustering analysis indicates 
that they are nonetheless underabundant relative to the expected 
number density of galaxies near QSOs. 

One of the observations that supports an intrinsic origin for 
many associated absorbers, is their solar or super-solar metallicity 
(see the Introduction). We do not find any evidence for such ele- 
vated metallicities in our sample of PDLAs. The same is true of 
the associated Mgll systems (some of which may be DLAs) stud- 
ied by vanden Berk et al. (2008). Nonetheless, as discussed above, 
we do find that PDLAs with high HI column densities have a mean 
metallicity that is higher than intervening systems by around a fac- 
tor of three. For example, [S/H]=-0.88 ± 0.24 for the proximate 
absorbers versus -1.41 ± 0.20 for those at AV > 10,000 km s~\ 
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QSO 


log N(HI) 


AVhi 


Velocity offset 


High 






(kms" 1 ) 


(kms^ 1 ) 


Ions 


J0142+0023 


20.38±0.05 


1772 


-100 


CIV, SilV 


Q0151+048 


20.34±0.02 


-1199 


-825 


NV, CIV 


J1131+6044 


20.50±0.15 


2424 


-800 


CIVOVI 


J1240+1455 


21.3±0.2 


102 


-125 


NV 


J 1604+3951 


21.75+.0.2 


-656 


-400 


CIV, SilV, NV 


J2321+1421 


20.70±0.05 


-1616 


-100 


CIV, SilV 



Table 10. High ionization species offset from low ion detections. No offset high ions are detected for J0140-0839 



The correlation between stellar mass and gas-phase metallicity that 
exists out to redshifts of at least 3 (Tremonti et al. 2004; Savaglio et 
al. 2005; Erb et al. 2006; Maiolino et al 2008) indicates that the PD- 
LAs might therefore also be relatively massive. QSOs themselves 
are highly biased tracers of mass and are often located in or near 
galaxy clusters or other overdensities, although they apparently es- 
chew the centres of clusters (Sanchez & Gonzalez-Serrano 1999; 
Barr et al. 2003; Sochting, Clowes & Campusano 2004; Kauffmann 
et al. 2004; Lietzen et al. 2009; Hutchings, Scholz and Bianchi 
2009). Since the typical stellar mass of a galaxy tends to increase 
with the galaxy density of its local environment (e.g. Baldry et al. 
2006), the combination of clustering and the mass-metallicity re- 
lation may explain the higher metallicities of PDLAs. Further evi- 
dence for this scenario comes from the discovery of a likely com- 
panion QSO to Jl 13 1+6044, and with the previously known com- 
panions of Q0151+048 and Q2343-BX415 (M0ller et al. 1998; 
Rix et al. 2007). Boris et al. (2007) find evidence that 3/4 of the 
QSO binaries in their sample are associated with rich clusters at 
z ~ 1. Other evidence for absorbers clustered around QSOs comes 
from the incidence of transverse absorbers. In their study of close 
projected pairs (in which it is argued that transverse absorbers suf- 
fer less photo-evaporation than line of sight absorbers) Hennawi et 
al. (2006) find a 50% covering fraction of sub-DLAs within 150 
kpc. 

It may seem suprising that absorbers with velocity offsets of 
up to 3000 km s -1 could be considered associated with the QSO 
environment when this velocity corresponds to a Hubble flow dis- 
tance of ~ 10 Mpc (proper) at z ~ 3. However, the velocity offsets 
are unlikely to indicate Hubble flow distances, as clearly demon- 
strated by the existence of absorbers with very negative velocities. 
Although we have additionally highlighted the uncertainty in the 
emission redshifts, these are unlikely to be incorrect by the > 1500 
km s -1 required to make all of our absorber offsets positive. We 
have discussed above some of the observations that have shown 
that QSOs often occupy overdense regions, a statement that might 
lead us to imagine cluster-like structures today with sizes of a few 
Mpc. However, Lietzen et al. (2009) have argued that QSOs can 
also trace out other overdense structures, such as filaments whose 
dimensions can be tens of Mpc. Haines et al. (2004) have found evi- 
dence for an extreme structure of galaxies clustered around quasars 
on scales of tens of Mpc at z ~ 1. 

The PDLAs appear to share some characteristics in common 
with absorbers associated with gamma-ray burst (GRB) hosts. The 
GRB DLAs have higher metallicities and are skewed to higher 
Af(Hi) than intervening QSO DLAs (e.g. Prochaska et al. 2007a). 
Higher [a/Fe] are also found in the GRB DLAs, similar to some 
of the PDLAs in our sample. However, the difference between the 
GRB DLAs and intervening systems is often attributed to the GRB 
sightlines preferentially intersecting the galaxy at smaller galacto- 



centric radii, rather than a fundamental difference in the populations 
probed by the GRBs and QSOs. 

Regardless of their origin, if the metallicities of (at least the 
high N(H i)) proximate DLAs are relatively high, as a population 
they present an interesting new selection technique for identify- 
ing the most metal-rich galaxies at high redshift. With metallici- 
ties ranging from 1/3 to 1/10 of the solar value, the high JV(Hi) 
PDLAs have metallicities similar to the 'metal strong DLAs' (MS- 
DLAs, Herbert-Fort et al. 2006) at z ~ 2 (Kaplan et al., in prepara- 
tion). Indeed, a number of the MSDLAs are also PDLAs (Kaplan et 
al., in preparation). Although this is still more metal-poor than the 
majority of Lyman Break Galaxies (LBG) studied at this redshift, 
most LBG abundances are limited to fairly massive galaxies (Erb et 
al. 2006). Gravitational lensing permits studies of fainter galaxies. 
Only a handful of such objects have been studied so far, but the re- 
sults indicate that galaxies with ~ luminosities at z ~ 2 to 3 have 
metallicities of Z ~ 1/2Z Q (e.g. Teplitz et al. 2000; Pettini et al. 
2002; Quider et al. 2009a,b). If the PDLAs follow a similar mass- 
metallicity relation, then they may be only slightly less massive 
than these lensed LBGs. High metallicity systems present a num- 
ber of interesting possibilities for studying the high redshift ISM. 
For example, Noterdaeme et al. (2008) have shown that molecu- 
lar hydrogen is highly dependent on metallicity. High metallicities 
also present the opportunity to detect and study the abundances of 
rarely detected elements such as Ge, B, CI and Co (Ellison, Ryan 
& Prochaska 2001; Prochaska et al. 2003c). 



9.4 Future Work 

Hennawi et al. (2006) presented a sample of close, projected QSO 
pairs where the background quasar exhibits optically thick Lyo- ab- 
sorption at the redshift of the foreground quasar. It is shown that 
50% of projected QSO pairs have an absorber with log A'(Hi) >19 
at the redshift of the foreground QSO when the separation is < 
150 kpc. This much higher incidence of transverse absorbers, rela- 
tive to proximate line of sight absorbers, is suggested by Hennawi 
& Prochaska (2007) to be due to differential photo-evaporation by 
anisotropic radiation of the foreground QSO. In the future, it will 
be interesting to compare the chemical abundances and ionization 
indicators of PDLAs to transverse DLAs at a given A^Hi) (e.g. 
Prochaska & Hennawi 2009) and to test fine structure distance es- 
timates against the measurable transverse separations. 

Despite their external provenance, the QSO's radiation can ap- 
parently affect PDLAs out to at least 2500 km s _1 from the QSO. 
With a larger sample, and accurate redshifts, it will also be possi- 
ble to look for trends with velocity. Ultimately, this may require the 
incorporation of other parameters such as QSO luminosity, includ- 
ing radiation from nearby companions. A campaign is currently 
underway to obtain IR spectra for a sample of QSOs with PDLAs 
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(including those presented in this paper) that will yield redshifts 
accurate to ~ 40 km s _I from the [OIII] A 5007 line. It is worth 
noting that only a few tens of z > 2 QSOs have redshifts [OIII] 
A 5007 redshifts (e.g. Scott et al. 2000). For objects in the Hub- 
ble flow, velocities correspond to distances. For positive values of 
AV we can not distinguish between these expansion velocities and 
peculiar motions (as expected if the galaxies and QSOs are in the 
same gravitational potential. However, negative AV can not be due 
to Hubble expansion, since it would imply the absorber is more dis- 
tant than the QSO, which is clearly impossible. The negative side 
of the velocity distribution of PDLAs therefore yields the 'true' 
velocity distribution of the PDLAs around the QSO, without con- 
tamination from the Hubble flow. Larger samples will also allow us 
to study trends of the combined effect of AV and QSO luminosity. 



10 CONCLUSIONS 

We have presented new high resolution echelle spectra for seven 
proximate damped Lyman alpha systems with AV < 3000 km s" 1 . 
The metal column densities derived from Voigt profile fits and the 
apparent optical depth method are complemented with abundances 
for a further 9 PDLAs taken from the literature. Our PDLA sam- 
ple is compared to the most complete sample of intervening DLAs 
currently available. Our principal conclusions are: 

(i) PDLAs exhibit a range of metallicities at a given redshift, 
ranging from ~ 1/3 to 1/1000 of the solar value (Sec tion [7, H and 
Figure [T4l. One of the PDLAs in our sample exhibits the lowest 
N(SiII)/N(HI) of any known DLA and has a value similiar to the 
intergalactic medium at this redshift. 

(ii) Based on this modest-sized sample, there is a general trend 
(with the exception of one fairly high velocity PDLA) for low 
metallicities (Z ~ l/50Z o ) in PDLAs with log N(Hi)< 21 and 
higher metallicities (Z ~ l/10Z o ) when the HI column density is 
higher (Figure [T5V At these high Hi column densities, the metal- 
licities of PDLAs are systematically higher than the intervening 
sample by a factor of around three. 

(iii) At least half of the PDLAs with A'(Hi) < 20.8 have sub- 
solar ratios of [S/Si] which can not be easily explained by known 
dust or nucleosynthetic trends (Section [6.3| and Figure [131 . We in- 
terpret the low values as resulting from ionization by a hard spec- 
trum. Sub-solar values of [S/Si] can be present even in PDLAs with 
AV > 2000 km s" 1 , and with no obvious trend with velocity. 

(iv) In addition to the dependence of metallicity on N(H i) and 
the sub-solar [S/Si] ratios, other indications of enhanced/hard ion- 
ization in the PDLAs which distinguish them from the intervening 
DLAs are: 1) a possibly higher fraction of NV absorbers (tenta- 
tive based on the small number statistics of our sample), 2) higher 
fractions of SilV/Sill at low A'(Hi), 3) similar velocity structure in 
SilV as Sill in 3/5 QSOs where the comparison can be made, 4) 
low ratios of [Ar/Si] (Section[6). 

(v) Most of the PDLAs (6/7) in our sample of seven have ad- 
ditional high ionization gas at large negative velocities of -100 
to -825 km s _I (Section 19.21 and Table 1 10b . The most extreme 
examples both have either a confirmed (Q0151+048) or tentative 
(11131+6044) close companion QSO. 

(vi) One of the PDLAs in our sample has very sub-solar [Fe/Zn] 
consistent with a large dust depletion fraction. However, in gen- 
eral, the [Fe/Zn] ratios of the PDLAs are consistent with the inter- 
vening population (Section[772]and Figure |T6t. The PDLA towards 
J 1604+ 3951 has very different depletion in its two main compo- 
nents. 



(vii) The range of alpha-to-iron ratios in the PDLAs is also con- 
sistent with the intervening DLAs, although there are a few notably 
high values (Section|73]and Figure [771. 

(viii) From analyses of the fine structure lines of Sill and CII, 
QSO-absorber distances at least a few tens of kpc are determined, 
with one case being constrained to have a separation of > 160 kpc 
(Section[8]l. 

We conclude that PDLA properties are generally consistent 
with an origin external to the QSO host (in contrast with, e.g. the 
narrow line associated systems, D'Odorico et al. 2004). However, 
the larger abundances (at least of high N(H i) PDLAs) imply that 
they may not be representative of the intervening sample. We sug- 
gest that the PDLAs may preferentially sample overdense environ- 
ments where biased galaxy formation has assembled more massive 
galaxies with higher metallicities. Indeed, at low redshift, some as- 
sociated absorbers have been identified with small clusters of galax- 
ies (e.g. Bergeron & Boisse 1986; Hamann et al. 1997a). If con- 
firmed with larger larger samples, this means that PDLAs could be 
used as probes of massive galaxies at high redshift. 

PDLAs have been widely excluded from most DLA surveys. 
The results presented here indicate that, depending on the science 
objective, this is a valid approach. We have argued that PDLAs 
may sample a rather special population of DLAs, possibly those 
in proto-clusters. Although this renders the PDLAs an interesting 
probe of high redshift galaxies, our results reveal some of the biases 
that could be introduced into statistical surveys that do not impose 
a AV limit in their DLA selection. A high priority for future work 
will be the improvement of emission redshift determinations. The 
results presented here imply that even at fairly large relative veloc- 
ities ionization may affect abundance determinations. Determining 
more accurate AV values in larger PDLA samples will provide an 
empirically motivated cut-off for studies of intervening DLAs. 
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